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Abstract 

A  dual  wavelength.  Q-switched  laser  cavity  has  been  successfully  designed  and 
assembled  to  evaluate  the  new  neodymium  (Nd)  based  materials.  Initial  characteri¬ 
zation  has  been  achieved  for  Nd:YAG  and  Nd:YLF.  The  results  indicate  that  for  a 
fixed  pulse  repetition  frequency,  the  delay  time  can  be  used  to  adjust  the  relative 
energy  between  the  1.06/im  and  1.32/im  pulses.  Any  deficiency  in  the  performance 
of  one  transition  can  be  made  up  by  performance  in  the  other  simply  by  changing 
the  relative  amount  of  time  the  population  inversion  is  allowed  to  build  up  for  each 
pulse. 

The  best  performance  was  obtained  using  a  13cm  cavity  length,  the  shortest 
possible  cavity  with  the  available  equipment.  The  gain  in  the  Nd:YLF  1.047/im 
cavity  was  so  high,  the  diffraction  efficiency  of  the  acousto-optic  Q-switch  was  insuf¬ 
ficient  to  hold  off  lasing  with  a  5%  ouput  coupler.  The  transmission  of  the  1.047/im 
output  coupler  used  for  dual  wavelength  operation  of  Nd:YLF  was  30%  whereas  the 
transmission  of  the  1.06/im  output  coupler  necessary  to  achieve  dual  wavlength  op¬ 
eration  in  Nd:YAG  was  5%.  At  a  pulse  repetition  frequency  for  each  wavelength 
of  2kHz,  and  with  820mW  of  absorbed  pump  power,  Nd:YLF  produced  pulse  ener¬ 
gies  of  40±2/iJ  and  15.3±0.6/iJ,  and  pulse  lengths  of  44±lns,  and  410±10ns  at  the 
wavelengths  of  1.047/im  and  1.32/im  respectively.  Nd:YAG  produced  pulse  energies 
of  33±l/i.J  and  27±l/tJ,  and  pulse  widths  of  32±0.4ns  and  183±lns  at  the  wave¬ 
lengths  of  1.06/im  and  1.34/im  respectively.  A  separate  set  of  dual  wavelength  data 
is  also  prented  for  Nd:YLF  at  5.0kHz. 
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I.  Introduction 

1.1  Background, 

1.1.1  Laser  radar  application.  The  basic  concept  of  a  laser  radar  involves 
sending  a  laser  pulse  out  to  a  target  and  sensing  the  reflected  energy.  There  are  a 
number  of  different  measurements  which  can  be  made  using  a  laser  radar  including 
the  intensity  image  of  a  scene  (which  is  the  equivalent  to  using  the  laser  as  a  flood¬ 
light),  the  range  from  the  sensor  to  the  target,  and  depending  on  the  type  of  laser 
system,  even  the  target’s  velocity.  The  application  that  this  work  is  directed  toward, 
involves  both  measuring  an  intensity  image  and  the  range  to  the  target.  Using  a 
single  pulse  to  build  the  entire  image  of  a  target  at  a  long  range  requires  a  large  laser 
which  can  produce  high  energy  pulses.  An  alternative  to  the  single  pulse  method 
uses  a  step-stare  technique  where  the  laser  system  focuses  down  on  a  small  portion 
of  the  target  and  measures  the  reflected  energy  and  range  on  just  that  small  portion. 
Building  the  image  is  done  by  scanning  over  the  entire  scene.  The  amount  of  energy 
required  for  this  last  technique  is  much  less  than  would  be  required  to  build  the  im¬ 
age  with  a  single  pulse.  The  pixel  registered  range  information  from  this  step-stare 
image  can  be  used  by  an  automated  system  to  distinguish  between  different  types 
of  targets  such  as  trees,  missile  launchers,  tanks,  and  trucks. 

One  laser  imaging  application  of  interest  to  the  Air  Force  is  detecting  targets 
which  are  either  camouflaged  or  buried  in  background,  with  a  sensor  system  mounted 
on  the  front  of  a  smart  munition.  This  type  of  munition  is  typically  designed  to 
be  autonomous  which  means  the  system’s  electronics  would  have  to  process  the 
data  from  the  sensor  package  in  order  to  find  and  track  the  target.  Automatically 


1-1 


processing  this  tremendous  amount  of  image  data  efficiently  requires  some  target 
feature  which  the  system  can  use  to  distinguish  the  target  from  its  surroundings.  One 
possible  feature  is  reflectivity  in  the  near  infra-red  (IR).  It  has  been  shown  that  the 
reflectivity  of  natural  objects  tends  to  decrease  as  the  wavelength  increases  whereas 
the  reflectivity  of  man  made  objects  remains  relatively  constant  or  increases  with 
wavelength  (1:95;  2:41).  Comparing  two  images  made  at  the  different  wavelengths, 
natural  objects  appear  relatively  darker  in  the  longer  wavelength  image,  whereas 
man  made  objects  appear  relatively  the  same  or  brighter. 

An  imaging  sensor  for  a  smart  munition  could  take  advantage  of  these  prop¬ 
erties  to  extract  man  made  objects  from  the  natural  background.  The  design  of 
the  laser  system  for  such  a  sensor  must  consider  several  factors.  The  laser  needs  to 
produce  two  wavelengths  in  order  to  accomplish  the  discrimination,  but  the  pulse 
repetition  rate  must  be  very  fast  to  achieve  the  necessary  image  frame  rate  for  both 
wavelengths.  The  image  frame  rate  needs  to  be  fast  enough  to  update  the  relative 
position  of  the  target  during  terminal  homing  so  the  munition  can  accurately  steer 
for  an  intercept.  Increasing  the  pulse  repetition  rate  too  far  causes  the  pulse  energy 
to  decrease  (3:623)  which  in  turn  decreases  the  maximum  range  of  the  sensor.  How¬ 
ever,  the  sensor  must  be  able  to  image  out  to  long  ranges  because  a  munition  has 
only  limited  maneuverability,  and  so  large  course  corrections  need  to  be  made  a  long 
time  before  target  impact. 

Increasing  the  pulse  repetition  rate  also  increases  the  pulse  length  of  the  laser, 
but  short  pulse  lengths  are  needed  to  make  accurate  range  measurements  for  target 
recognition  purposes.  A  range  precision  on  the  order  of  one  foot  is  needed  to  make 
an  accurate  measurement  of  the  shape  of  a  target.  This  range  precision,  which  is  a 
measurement  of  the  difference  in  range  between  two  adjacent  pulses,  requires  a  laser 
pulse  length  on  the  order  of  10  nanoseconds. 

The  design  of  a  sensor  for  the  front  of  a  smart  munition  must  be  very  compact 
and  efficient.  To  meet  the  stringent  size  and  efficiency  constraints,  the  ideal  laser 
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design  would  use  a  single  laser  to  produce  both  wavelengths.  The  difficulty  is  finding 
a  suitable  laser  crystal  which  can  produce  two  wavelengths. 

The  maturity  of  the  technology  plays  an  important  role  when  selecting  a  laser 
crystal.  At  present,  the  most  mature  laser  technology  which  can  be  packaged  into 
the  size  requirements  necessary  for  an  autonomous  sensor  are  crystals  doped  with 
the  rare  earth  ion,  neodymium  (Nd).  These  crystals  can  be  diode  pumped  and  lasers 
have  been  constructed  with  this  technology  which  are  under  three  inches  long  and 
produce  average  powers  of  three  Watts  (5:89). 

1.1.2  Review  of  laser  cavity  designs.  As  pointed  out  earlier,  short  laser 
pulse  lengths  are  needed  to  make  accurate  three-dimensional  measurements  of  a  tar¬ 
get.  Considerable  progress  has  been  made  in  recent  years  to  develop  lasers  based  on 
neodymium  which  can  produce  short  pulses.  The  applications  currently  driving  this 
development  are  laser  communications  and  laser  imaging  which  need  pulse  lengths 
on  the  order  of  a  few  nanoseconds  or  less.  Pulse  lengths  this  short  have  been  pro¬ 
duced  in  several  neodymium  based  materials  pumped  with  laser  diodes  with  much 
success  (3:622;  20:1131;  21:6616;  22:885). 

One  particular  group  of  workers  from  Lightwave  Electronics,  Inc.  has  made 
significant  advances  in  producing  short  pulses,  and  the  account  of  their  progress  is 
detailed  in  a  series  of  articles  which  began  in  1990.  In  the  first  article,  they  outline 
a  patented  laser  cavity  design  which  has  a  very  high  gain  and  a  very  short  photon 
lifetime(3:622)  -  two  characteristics  essential  for  producing  short  pulses  (23:2044). 
This  group  achieved  a  short  photon  lifetime  by  having  a  very  short  cavity  which  is  less 
than  an  inch  long.  What  makes  this  group’s  achievements  all  the  more  interesting 
is  that  they  use  an  acousto-optic  modulator  to  Q-switch  their  laser. 

An  acousto-optic  device  deflects  light  by  sending  an  acoustic  wave  through  a 
quartz  or  glass  crystal.  Acousto-optic  devices  are  less  expensive  than  electro-optic 
devices,  which  are  the  only  other  alternative  for  applications  needing  controllable  Q- 
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switching.  Besides  cost,  electro-optic  devices  also  require  high  voltages  to  operate. 
The  disadvantage  of  acousto-optic  Q-switches  is  that  they  take  longer  to  switch  which 
often  results  in  longer  pulses.  The  pulse  lengths  the  authors  were  able  to  achieve 
were  on  the  order  of  5-8ns  at  a  wavelength  of  1.06/im  and  25-40ns  at  a  wavelength 
of  1.32/im. 

The  second  article  from  this  group,  published  in  1991  (22:885),  outlines  how 
they  achieved  pulse  lengths  of  2.7ns  at  repetition  rates  of  up  to  10kHz  in  Nd:YLF 
(neodymium  yttrium  lithium  fluoride)  with  a  laser  cavity  about  a  quarter  inch  long. 
These  pulse  lengths  are  on  the  same  order  needed  to  make  range  measurements  with 
a  range  resolution  on  the  order  of  one  foot. 

The  group’s  last  article  was  published  in  1993  (21:6616)  and  describes  work 
with  the  same  cavity  design  as  the  one  used  in  1990  but  with  an  additional  laser  crys¬ 
tal  material,  Nd:YV04  (neodymium  yttrium  orthovanadate).  The  technical  achieve¬ 
ment  they  highlight  in  this  article  as  being  crucial  to  their  success  was  an  acousto¬ 
optic  Q-switch  that  had  a  switching  time  of  18ns.  Using  their  design,  they  were  able 
to  achieve  0.6ns  pulses  at  a  repetition  rate  of  1kHz. 

Although  all  of  these  short-pulse  lasers  produced  relatively  low  powers,  a  dif¬ 
ferent  group,  Bear  et  al.  (20:1131),  has  developed  a  method  for  producing  much 
higher  powers  using  an  alternate  cavity  design.  They  report  the  highest  efficiency 
ever  in  a  diode-pumped  solid-state  laser  operating  in  the  TEMqo  mode.  Using  the 
materials  Nd:YAG  (neodymium  yttrium  aluminum  garnet)  and  Nd:YLF,  they  were 
able  to  achieve  pulse  lengths  of  5ns  at  pulse  repetition  frequencies  of  1kHz.  Recently, 
working  at  Spectra-Physics,  Bear  has  produced  a  Nd:YV04  laser  which  produces  5-6 
Watts  of  average  power.  The  pulse  repetition  frequency  for  this  laser  is  20kHz  with 
pulse  lengths  on  the  order  of  7-8ns  (9).  In  order  to  apply  their  advances  to  the 
laser  of  interest,  a  method  for  producing  two  wavelengths  with  their  cavity  design  is 
required. 
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All  of  the  lasers  discussed  above  which  produce  short  pulse  lengths  are  single 
wavelength  devices.  In  order  to  change  wavelengths,  it  is  necessary  to  also  change 
the  optics.  There  are,  however,  several  cavities  already  developed  which  can  pro¬ 
duce  two  wavelengths  from  the  same  crystal.  Simultaneous  lasing  at  both  1.06/zm 
and  1.34/im  wavelengths  has  been  reported  in  several  neodymium  doped  crystals  in¬ 
cluding  Nd:YAG,  Nd:YLF,  Nd:BEL,  NdYAP  (24:2315;  25:746;  26:2588;  27:444)  and 
most  recently  Nd:SFAP  (32).  The  work  in  this  area  has  been  driven  by  an  application 
which  needs  simultaneous  pulses  at  both  wavelengths  for  sum  harmonic  frequency 
generation.  An  imaging  laser  radar  would  benefit  from  simultaneous  pulses  because 
any  tracking  error  would  be  essentially  eliminated  allowing  the  reflectivities  for  a 
single  target  pixel  to  be  compared.  Pulses  separated  in  time  may  require  extra  pro¬ 
cessing  to  ensure  that  the  intensities  being  compared  correspond  to  the  same  point 
on  the  target.  If  two  pulses  of  different  wavelengths  were  sent  out  to  a  target  simul¬ 
taneously,  a  method  of  separating  them  in  the  receiver  would  be  necessary  to  process 
the  information  they  contain.  The  work  in  this  area  relevant  to  the  laser  of  interest 
is  the  “Y”-cavity  laser  cavity  design  (see  Figure  1.1)  developed  by  Tomaschke  et  al. 
(33:252)  for  producing  two  wavelengths  from  a  single  crystal.  A  dichroic  beamsplitter 
is  used  to  separate  the  two  wavelengths  by  reflecting  one  wavelength  and  transmit¬ 
ting  the  other  into  different  legs,  each  containing  an  acousto-optic  Q-switch.  The 
authors  develop  a  model  for  this  system  but  it  is  only  applicable  for  simultaneous 
wavelength  pulses,  not  alternating  pulses. 

1.2  Problem 

There  is  an  ongoing  quest  for  the  ultimate  host  material  for  neodymium  which 
will  yield  a  combination  of  high  pumping  efficiency,  high  gain,  thermal  stability,  and 
good  optical  quality  (12:387).  This  quest  has  generated  a  number  of  attractive  crys¬ 
tals.  The  difficulty  arises  in  choosing  the  best  laser  crystal  for  this  dual  wavelength 
laser  radar  application.  A  numerical  analysis  to  choose  the  best  crystal  would  have 
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Figure  1.1  Tomaschke’s  laser  cavity  configuration  for  simultaneous  wavelength  out¬ 
put  (pulsed). 

to  rely  on  the  published  data  for  the  different  material  parameters  of  each  of  the 
available  laser  crystals.  Unfortunately,  either  these  reported  values  vary  tremen¬ 
dously,  or  the  values  have  not  been  published.  The  only  way  to  draw  meaningful 
conclusions  about  the  relative  performance  of  these  available  crystals  is  to  measure 
their  performance  in  an  identical  laser  setup. 


1.3  Scope 

Equations  describing  continuous  (CW)  and  Q-switched  laser  performance  are 
presented  and  a  general  figure  of  merit  is  discussed  for  choosing  the  best  laser  ma¬ 
terial  based  on  published  data.  A  dual-wavelength  Q-switched  laser  is  designed 
and  constructed  which  has  a  “Y”  shaped  cavity.  The  performance  of  the  system  is 
characterized  for  two  different  laser  crystals,  Nd:YAG  and  Nd:YLF.  Both  CW  and 
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Q-switched  measurements  were  made.  The  CW  measurements  include  laser  pump 
threshold  and  slope  efficiency.  The  Q-switched  measurements  were  made  for  both 
single  and  dual  wavelength  operation,  and  include  the  laser  pulse  energy  and  pulse 
length  as  a  function  of  three  variables:  1)  pump  power,  2)  pulse  repetition  frequency, 
and  3)  the  delay  time  between  the  two  pulses.  A  numerical  model  is  fitted  to  the 
data  from  the  1.06/xm  results. 

1.4  Approach 

In  Chapter  II,  a  description  of  the  lasing  mechanism  in  Nd:YAG  provides  the 
foundation  for  presenting  the  equations  describing  both  CW  and  Q-switched  laser 
performance.  Using  these  concepts,  a  figure  of  merit  is  discussed  for  choosing  the 
most  promising  material  from  data  in  the  published  literature. 

Laser  crystals  made  from  Nd:YAG,  Nd:YLF,  Nd:YV04,  Nd:GdV04  (neodymium 
gadolinium  orthovanadate),  and  Nd:SVAP  (neodymium  strontium  fluorovanadate) 
were  purchased  for  the  system  but  only  Nd:YAG,  Nd:YLF,  and  Nd:SVAP  arrived 
in  time  for  evaluation.  In  Chapter  III,  measurements  of  the  polarized  absorption 
spectra  of  each  material  in  the  range  from  790  nm  out  to  815  nm  are  presented,  and 
the  measurements  indicate  the  best  wavelength  for  pumping  in  the  4Hg/2  to  4Fs/2 
band. 

A  “Y”  shaped  cavity  was  constructed  which  allows  Q-switched  operation  at 
alternating  wavelengths  (see  Figure  1.1).  Initial  laser  measurements  include  finding 
the  CW  laser  threshold  for  both  the  nominal  1.06/im  and  1.32/xm  transitions  in  each 
sample.  The  change  in  threshold  caused  by  inserting  each  component  was  used  to 
calculate  an  estimate  of  the  intra-cavity  loss  caused  by  each  component  at  the  lasing 
wavelength. 

Q-switched  laser  performance  measurements  include  the  pulse  energy  and  pulse 
length  (FWHM)  as  functions  of  pulse  repetition  frequency.  These  measurements 
were  made  for  two  different  cavity  lengths.  After  comparing  these  results  against 
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theoretical  predictions  from  the  model,  both  the  measurement  and  model  results  are 
used  to  draw  conclusions  about  the  best  material  and  cavity  configuration  for  the 
laser  imager  application. 


1-8 


II.  Laser  Theory 


2.1  Population  inversion  and  gain  in  neodymium  (Nd)  based  crystals 

Figure  2.1  shows  the  energy  band  diagram  for  Nd  in  YAG  (yttrium  aluminum 
garnet)  (6:49).  Although  the  details  of  this  diagram  change  when  Nd  is  placed  in 
different  host  crystals,  the  basic  electronic  transitions  which  produce  lasing  remain 
the  same.  A  pump  photon  from  a  laser  diode  with  a  wavelength  around  808nm 
can  cause  an  electron  in  the  4Ig/2  ground  state  to  transition  to  the  band  of  allowed 
states  between  the  4Hg/2  and  the  4F5//2  states.  The  width  of  this  band  determines  the 
absorption  bandwidth,  the  range  of  pump  wavelengths  that  can  be  used  to  pump 
the  laser.  High  power  laser  diode  arrays  experience  relatively  large  differences  in 
the  wavelength  from  the  many  elements  which  make  up  the  arrays  and  so  wider 
absorption  bandwidths  allow  increased  absorption  of  the  pump  energy  (8:571;  10:999; 
15:298).  The  most  likely  decay  route  from  this  band  of  allowed  states  between 
the  4Hg/2  and  the  4F5/2  levels  is  into  the  4F3/2  state  via  a  non-radiative  relaxation 
process  where  the  energy  is  given  up  into  lattice  vibrations.  The  4F3/2  state  is  the 
upper  energy  level  for  both  the  transitions  which  produce  the  1.06/xm  and  1.32/im 
wavelengths.  The  average  amount  of  time  that  electrons  spend  in  this  state  is  referred 
to  as  the  upper  state  lifetime,  tj. 

If  an  electron  in  the  upper  laser  energy  level  transitions  down  to  the  4In/2 
state,  a  photon  with  a  wavelength  of  1.06  pm  will  be  emitted.  If  this  emitted  photon 
interacts  with  another  neodymium  ion  which  has  an  electron  in  the  upper  laser 
transition  level,  it  could  induce  that  electron  to  transition  to  the  4In/2  state,  causing 
another  1.06  pm  photon  to  be  emitted.  This  process  is  called  stimulated  emission. 
A  measure  of  the  probability  that  this  stimulated  emission  process  will  take  place  is 
the  stimulated  emission  cross  section,  ctse- 

A  related  process  is  called  stimulated  absorption.  This  process  also  has  an 
associated  cross  section,  a  Ah-  Absorption  can  occur  when  a  photon  encounters  an 


2-1 


Pump 

band 


^9/2 

- 4F5/; 

_ **3/5 

Pump 
10 

Energy 
(lO^cm'1) 

_ *13/2 

4 

- Ill/2 

4 

I9/2 

Figure  2.1  Energy  band  diagram  for  neodymium  in  YAG. 

ion  with  an  electron  in  the  lower  4In/2  energy  level.  The  ion  aborbs  the  energy  in 
the  photon  causing  a  transition  to  the  upper  energy  level. 

Population  inversion  occurs  when  the  number  of  electrons  in  the  4F3/2  level 
exceeds  the  number  in  the  lower  energy  level.  This  situation  implies  that  a  photon 
traversing  the  laser  crystal  is  more  likely  to  encounter  an  electron  in  the  upper  energy 
level  than  it  is  to  encounter  an  electron  in  the  lower  laser  energy  level.  Stimulated 
emission  is  then  more  probable  in  this  case  than  absorption,  and  so  for  a  single  pass 
through  the  material,  a  spontaneously  emitted  photon  will  see  a  net  gain.  At  room 
temperature  in  thermal  equlibrium,  most  of  the  electrons  reside  in  the  4I9/2  state. 
The  situation  where  more  electrons  reside  in  the  upper  energy  state  is  not  a  natural 
occurance  and  thus  a  pumping  mechanism  is  required  to  achieve  it. 
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A  laser  crystal  being  pumped  with  a  constant  source  has  a  gain  per  unit  length 
in  the  gain  medium  given  by(35:183) 


7o  —  &SE 


(2.1) 


where  use  is  the  stimulated  emission  cross  section  for  the  lasing  transition,  jV2  is 
the  population  density  of  the  upper  energy  state,  N\  is  the  population  density  of  the 
lower  energy  state  for  the  lasing  transition,  and  g2  and  gi  are  the  degeneracies  of  the 
upper  and  lower  laser  energy  levels  respectively.  The  gain  coefficient  is  positive  for 
the  case  where  iV2  >  ^N\.  When  iV2  <  the  gain  coefficient  is  negative,  and  a 
spontaneously  emitted  photon  experiences  a  net  loss. 


In  the  absence  of  any  feedback,  the  population  density  in  state  2  builds  ac¬ 
cording  to  (35:191) 


JV2(t)  =  R2tj 


(2.2) 


where  R2  is  the  pumping  rate  per  unit  volume.  The  expression  for  f?2  is  given  by 


r2  = 


VAbsPpump 

hUpVoo 


(2.3) 


where  Ppump  is  the  pump  power  in  Watts,  h  is  Plank’s  constant,  v  is  the  frequency  of 
the  pump  radiation,  rjAbs  is  the  probability  that  the  incident  photons  are  converted 
to  useful  upper  level  population,  and  Voo  is  the  pumped  volume  in  the  gain  medium. 
The  value  for  r/Abs  in  this  expression  includes  the  fraction  of  incident  pump  photons 
absorbed  by  the  gain  medium  multiplied  by  the  fraction  that  actually  relax  into  the 
upper  energy  level  for  the  lasing  transition.  The  spatial  profile  of  the  pump  beam  is 
modeled  as  a  Gaussian  and  so  the  expression  for  Voo  is  given  by  (35:118) 

7ru>2 

^00  =  —l,  (2.4) 
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where  lo  is  the  spot  size  of  the  laser  beam  inside  the  laser  crystal. 

Assuming  both  steady  state  conditions  =  0,  and  that  the  lower  energy  state 
decays  infinitely  fast  such  that  Nx  =0,  the  population  inversion  is  just  equal  to  the 
population  in  state  2.  The  lifetime  of  the  lower  energy  level  for  the  4F3/2  -  4In/2 
transition  has  been  measured  in  Nd:YAG  and  Nd:YLF(7:1987).  Nd:YAG  has  a  lower 
level  lifetime  of  0.225ns  and  Nd:YLF  has  a  lifetime  of  20ns.  Given  this  data,  the 
approximation  that  Ni  =0  is  valid  for  Nd:YAG  for  pulse  lengths  over  0.225ns  and 
for  Nd:YLF  for  pulse  lengths  over  20ns.  The  shortest  pulse  1.06/im  lengths  achieved 
in  this  work  were  35ns  long  (see  Chapter  IV)  and  so  the  assumption  that  Nj  =0  is 
valid  for  this  particular  set  up.  Using  these  assumptions,  substituting  Eq  (2.3)  and 
Eq  (2.2)  into  the  expression  for  the  small  signal  gain  in  Eq  (2.1)  yields 

Ppump  VAbs  /n  r\ 

^  =  °SETI1^VM  (2'6) 

From  the  dependence  of  the  gain  on  the  upper  state  lifetime  illustrated  in 
Eq  (2.5),  it  can  be  seen  that  materials  with  longer  upper  state  lifetimes  are  better 
at  storing  the  pump  energy.  When  the  pulse  repetition  rate  is  below  T,  materials 
with  longer  lifetimes  can  achieve  higher  gain,  given  a  constant  pump  rate  (3:623). 

Eq  (2.5)  also  shows  the  dependence  of  the  gain  on  the  stimulated  emission 
cross  section,  a  larger  cross  section  means  less  feedback  is  needed  to  cause  lasing  to 

occur  and  it  also  means  that  shorter  pulses  are  possible.  The  value  for  the  stimulated 

emission  cross  section  applies  to  a  single  transition  in  a  particular  material.  Some  of 
the  published  values  for  this  cross  section  in  Nd:YAG  at  1.064/mi  vary  up  to  ±  30%. 

The  expression  for  the  gain  also  illustrates  the  dependence  on  the  size  of  the 
pumped  volume  in  the  laser  crystal.  Focusing  the  pump  down  into  a  smaller  spot 
size  increases  the  gain. 


2-4 


Figure  2.2  Simple  end-pumped  cavity  configuration. 

2.2  CW  laser  theory 

If  a  population  inversion  is  created  in  a  laser  crystal  and  feedback  is  provided 
to  the  stimulated  emission  process,  the  number  of  photons  produced  rises  rapidly. 
This  feedback  can  be  created  by  placing  the  laser  crystal  inside  a  Fabry  Perot  cavity. 
Figure  2.2  shows  some  of  the  basic  parameters  for  a  simple  laser  cavity.  In  order 
for  laser  oscillation  to  occur,  the  round  trip  gain  must  be  greater  than  the  round 
trip  losses.  This  requirement  means  that  the  gain  coefficient  must  be  greater  than  a 
threshold  given  by  (35:184) 

lth  =  2Tln  (r  W  r)  (2,6) 

Ztg  \lllCompL2J 

where  lg  is  the  length  of  the  gain  medium,  Ti  and  T2  represent  the  reflectivity  of  the 
cavity  end  mirrors,  and  Tcomp  represents  the  single  pass  transmission  of  the  various 
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optical  components  inside  the  laser  cavity  including  the  laser  crystal  itself.  Combin¬ 
ing  Eq  (2.1)  for  the  gain  coefficient  and  Eq  2.6  for  the  threshold  gain  coefficient  yields 
an  expression  for  the  threshold  population  inversion  nth  which  is  given  by  (35:257) 


nth 


Mg^SE 


In 


,riT£ompr2< 


Voo 


(2.7) 


After  the  population  inversion  exceeds  the  threshold,  the  general  equation 
which  describes  the  output  characteristics  of  a  CW  laser  is  given  by  (6:96) 


P out  —  °S  {Ppump  Prhreshold) 


(2.8) 


where  Ppump  is  the  incident  pump  power,  ^Threshold  is  the  threshold  pump  power,  and 
as  is  the  slope  efficiency  of  the  output  power  vs  the  input  power.  The  expression  for 
the  slope  efficiency  is  given  as  (6:100). 


as  = 


(l  T 2)Vefficiency 


1  + 


Comp 


[i  - 


(2.9) 


where  T i  is  the  power  reflectivity  of  the  cavity  high  reflector,  T2  is  the  power  reflec¬ 
tivity  of  the  output  coupler,  Tcomp  is  the  one  way  transmission  of  the  components 
within  the  cavity  including  the  laser  crystal  itself,  T]efficiency  is  the  pumping  efficiency, 
and  Ax,  is  the  wavlength  of  the  lasing  transition.  The  expression  for  the  pumping 
efficiency  is  given  by 

Ap  .  . 

V efficiency  =  VBVAbs~{  (2.10) 

Ajr 

where  tjb  is  the  beam  overlap  efficiency  which  is  inserted  to  include  only  the  gain 
which  is  within  the  mode  volume,  and  Xp  is  the  wavelength  of  the  pump  source.  The 
laser  cavity  shown  in  Figure  2.2  shows  an  end-pumped  scheme.  The  size  of  the  pump 
beam  relative  to  the  size  of  the  oscillating  laser  mode  effects  the  pumping  efficiency 
because  the  portion  of  the  gain  medium  outside  the  laser  mode  volume  does  not  need 
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to  be  pumped.  The  optimum  design  condition  is  when  the  pump  beam  size  and  the 
laser  beam  size  match  inside  the  laser  crystal.  If  Gaussian  pump  and  laser  profiles 
are  assumed,  the  beam  overlap  efficiency  can  be  calculated  from  (6:92) 


Vb 


2co2 

coj  +u>2 


for  LOp  >  co 


(2.11) 


and 

tjb  =  1  for  cop  <  co  (2.12) 

where  co  is  the  spot  size  (radial  distance  from  the  center  of  the  beam  to  the  ^ 
intensity  point)  of  the  laser  mode  in  the  gain  medium  and  u )p  is  the  spot  size  of  the 
pump  beam  in  the  gain  medium.  Eq  (2.12)  indicates  that  the  overlap  efficiency  is 
100%  when  cop  <  co. 

The  expression  for  the  amount  of  pump  power  to  reach  threshold,  Prhreshoid, 
is  given  by  (6:100) 


^Threshold  — 


-boo huL  In  (^TiT^p) 

&  S  ET f  ^efficiency  lg 


(2.13) 


where  Voo  is  the  laser  mode  volume  in  the  crystal,  h  is  Planck’s  constant,  uL  is  the 
optical  frequency  of  the  laser  radiation,  r2  is  the  reflectivity  of  the  output  coupler,  Pi 
is  the  reflectivity  of  the  high  reflector,  gse  is  the  stimulated  emission  cross  section, 
Tf  is  the  fluorescence  lifetime,  lg  is  the  length  of  the  gain  medium,  and  rjefficiency  is 
the  probability  that  an  incoming  photon  is  converted  to  useful  population  inversion. 

By  measuring  the  threshold  pump  power  with  different  output  coupler  reflec¬ 
tivities,  an  estimate  of  the  round  trip  cavity  loss  and  the  pumping  efficiency  can  be 
estimated  by  plotting  the  natural  log  of  the  output  coupler  reflectivity  against  the 
threshold  pump  power.  The  equation  describing  this  plot  is  (6:101) 


-  ln(r2)  =  2 
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Figure  2.3  Simple  Q-switched  cavity  configuration. 


where  L  is  the  total  loss  in  the  cavity  other  than  the  output  coupler 


I  _  p  rp2 

1  1  I-1- camp 


(2.15) 


2.3  Q-switched  laser  theory 

Q-switching  a  laser  cavity  is  the  process  of  increasing  the  cavity  loss  to  the 
point  where  no  lasing  can  take  place.  Under  continuous  pumping,  this  causes  a 
large  population  inversion  to  be  created.  Then  the  loss  in  the  cavity  is  reduced  very 
quickly,  causing  feedback  to  the  gain  medium  to  increase.  The  stimulated  emission 
process  quickly  escalates,  resulting  in  a  large  number  of  photons  emitted  from  the 
laser  in  a  short  pulse,  all  with  roughly  the  same  wavelength. 

The  diagram  for  a  simple  Q-switched  cavity  configuration  is  shown  in  Fig¬ 
ure  2.3.  For  this  analysis  the  pump  is  treated  as  a  CW  optical  beam  from  another 
laser.  The  Q-switch  element  provides  additional  loss  in  the  cavity,  this  additional 
loss  causes  the  threshold  population  level  to  increase  according  to  Eq  (2.7).  If  this 
increased  threshold  level  is  high  enough,  the  actual  population  inversion  remains  be¬ 
low  this  threshold  and  the  cavity  does  not  lase.  Under  this  situation  of  CW  pumping 
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and  high  Q-switch  loss,  the  population  builds  up  according  to  Eq  (2.2).  After  the 
Q-switch  returns  to  low  loss,  the  photon  number  Nphot  and  the  population  inversion 
n  change  according  to  the  two  coupled  nonlinear  differential  equations  (35:253) 


dn  n 

~  ~~ Nphot  (2-17) 

dl  Tlth 

where  T  is  time  normalized  to  units  of  photon  lifetime,  and  nth  is  the  threshold 
population  inversion.  The  expression  for  the  photon  lifetime  is  given  by 


Tphot 


1  —  VlT2T^omp 


(2.18) 


where  rrt  is  the  round  trip  time  for  a  photon  in  the  cavity. 

The  two  differential  equations  which  describe  the  change  in  the  number  of 
photons  circulating  in  the  cavity  (Eq  (2.16)),  and  the  change  in  the  total  population 
inversion  (Eq  (2.17)),  assume  that  the  effect  of  the  pumping  mechanism  is  negligible 
over  the  duration  of  the  pulse  and  that  the  Q-switch  has  essentially  a  zero  turn-on 
time.  These  equations  were  implemented  in  a  software  model  (see  Appendix  B) 
and  Figure  2.4  on  page  2-10  shows  an  example  of  the  population  inversion  and  the 
photon  number  as  functions  of  time.  The  threshold  population  level  in  the  upper 
plot  changes  because  the  Q-switch  is  modeled  as  having  a  non-zero  switching  speed, 
and  the  additional  loss  caused  by  the  Q-switch  causes  a  change  in  the  threshold. 

The  initial  photon  number  and  initial  population  level  used  to  initiate  these 
equations  have  a  significant  effect  on  the  results.  The  photons  which  exsist  in  the 
cavity  prior  to  lasing  that  have  wavelengths  corresponding  to  the  lasing  transition, 
are  a  result  of  spontaneous  emission.  The  initial  number  of  these  photons  has  an 
indirect  impact  on  the  amount  of  energy  calculated  in  the  pulse  and  the  pulse  width 
by  effecting  on  how  fast  the  pulse  builds  up  after  the  threshold  population  level  drops 


Time  (ns) 

Figure  2.4  Population  inversion,  threshold  population  level,  and  photon  number  as 
a  function  of  time. 

below  the  initial  population  inversion  level.  When  the  Q-switch  is  modeled  as  taking 
a  finite  time  to  switch  from  high  loss  to  low  loss,  as  it  is  here,  a  large  initial  photon 
number  causes  the  pulse  to  build  up  earlier  while  the  Q-switch  is  still  switching  from 
high  to  low  loss.  This  overlap  of  the  Q-switch  transition  with  the  emerging  pulse 
causes  a  decrease  in  the  amount  of  energy  in  the  pulse  and  an  increase  in  the  pulse 
length.  A  lower  initial  photon  number  causes  the  pulse  to  build  up  later  when  the 
Q-switch  is  almost  fully  transmissive  and  has  no  significant  effect  on  the  pulse. 
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The  total  energy  in  the  output  pulse  shown  in  the  bottom  portion  of  Figure  2.4 
is  calculated  from  the  change  in  the  inversion  level  which  is  shown  in  the  top  portion 
of  Figure  2.4.  The  expression  for  the  laser  output  energy  is  given  by  (23:2042) 


Wont  =  (n*  -  nf)  huL — —  (2.19) 

Q-total 

where  n*  and  rif  are  the  initial  and  final  population  inversions  respectively,  aext  is 
the  portion  of  loss  which  is  extracted  from  the  cavity  as  useful  output,  and  atotai  is 
the  total  loss  in  a  round  trip  of  the  cavity.  The  expression  for  atotai  is  given  by 


atotai 


In 


,rir2T c0mp, 


(2.20) 


where  again,  Fi  is  the  power  reflectivity  of  the  high  reflector,  P2  is  the  power  reflec¬ 
tivity  of  the  output  coupler,  and  Tcomp  is  the  total  one  way  transmission  of  all  of 
the  various  components  in  the  cavity.  The  expression  for  aext 


&ext  —  1^ 


(2.21) 


A  relationship  exists  beween  the  initial  population  level  at  the  start  of  the 
pulse,  the  threshold  level,  and  the  fraction  of  energy  which  is  extracted  from  the 
population  inversion.  The  initial,  final,  and  threshold  inversion  levels  are  related 
through  the  expression  (35:255) 


(2.22) 


The  fraction  of  the  initial  inversion  which  is  converted  to  photons  is  typically  defined 
as  an  extraction  efficiency 


Vxtrn  — 


rii  -  rif 

Ui 


(2.23) 
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The  relationship  of  this  value  to  the  inversion  ratio,  the  ratio  of  initial  population 
inversion  n*  to  threshold  inversion  nth  is  shown  in  Figure  2.5.  From  this  relationship, 


Figure  2.5  Extraction  efficiency  as  a  function  of  the  inversion  ratio  ^LL. 

it  becomes  evident  that  the  percentage  of  the  energy  extracted  from  the  initial  pop¬ 
ulation  inversion  is  a  strong  function  of  the  inversion  ratio.  For  a  Q-switched  pulse, 
the  initial  population  level  is  very  important.  According  to  Eq  (2.2),  it  depends  on 
the  pump  rate  per  unit  volume,  the  fluorescence  lifetime,  and  the  amount  of  time  the 
population  is  allowed  to  build  up.  In  a  repetitively  Q-switched  laser,  another  factor 
to  consider  is  the  residual  population  inversion  from  the  previous  pulse.  This  leftover 
population,  depending  on  how  efficiently  the  last  pulse  swept  out  the  gain  region, 
can  be  significant.  In  order  to  properly  model  the  laser’s  behavior,  this  number  must 
be  taken  into  account  by  calculating  an  equivalent  pumping  time.  The  equivalent 
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pumping  time  is  the  amount  of  time  that  it  would  take  to  reach  this  final  population 
level  at  the  given  pumping  rate  and  it  can  be  calculated  by  using  the  final  population 
inversion  n/  from  the  last  pulse  in  Eq  (2.2)  and  solving  for  time 

nf  \ 

R*Tf) 

The  initial  population  inversion  is  then  calculated  using  the  sum  of  this  equiva¬ 
lent  pumping  time  with  pU^sg'r1ep  rate  in  for  t  in  Eq  (2.2).  A  model  of  the  Q-switched 
laser  performance  was  developed  using  the  equations  in  sections  2.1  -  2.3.  A  version 
of  this  model,  written  in  Matlab,  is  presented  in  Appendix  B. 

2.4  General  figure  of  merit  for  dual  wavelength  laser  performance 

The  performance  of  a  particular  laser  system  depends  on  a  large  number  of 
variables.  Examining  the  expression  for  the  threshold  pump  power  in  Eq  (2.13)  on 
page  2-7,  there  are  a  large  number  of  factors  which  are  based  on  the  particular  cavity. 
Lower  threshold  values  can  be  achieved  by  reducing  both  the  size  of  the  beam  in  the 
crystal  and  the  cavity  losses.  For  the  sake  of  comparing  different  materials,  these 
will  be  assumed  to  be  the  same  for  all. 

The  significant  material  dependent  parameters  in  Eq  (2.13)  are  the  absorp¬ 
tion  efficiency,  stimulated  emission  cross  section,  and  fluorescence  lifetime.  The 
absorption  efficiency  is  a  function  of  the  material  and  the  pumping  mechanism.  A 
Ti:Sapphire  laser  has  a  very  narrow  line  width  which  allows  this  source  to  be  a  very 
efficient  pump  for  most  materials  even  when  the  absorption  line  width  in  the  laser 
material  is  very  narrow.  High  power  laser  diode  arrays  used  for  laser  diode  pumping 
exhibit  some  wavelength  variation  among  the  individual  elements  (18:560).  In  order 
to  use  the  pump  power  most  efficiently,  the  absorption  line  width  in  the  material 
should  be  broader  than  the  wavelength  content  from  the  laser  diode  arrays.  Ab¬ 
sorption  efficiency  is  a  complex  problem  dependent  on  the  laser  material,  pumping 


1 'equiv 


=  —Tf  In  1 
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mechanism,  and  laser  geometry.  Since  the  total  absorptance  can  be  made  equal  in 
each  of  the  materials  by  changing  the  length  of  the  gain  region,  it  will  be  set  aside 
for  the  sake  of  comparison. 

The  remaining  factor  is  the  product  of  the  stimulated  emission  and  the  upper 
state  lifetime,  crsETf.  The  values  of  these  two  parameters  depends  on  the  interaction 
of  the  neodymium  ion  with  the  surrounding  crystal  lattice.  The  larger  this  product 
is,  the  lower  the  threshold,  and  so  this  product  provides  a  general  figure  of  merit 
for  comparing  different  materials  (12:159).  The  performance  of  a  laser  operating  at 
two  wavelengths  depends  on  this  product  from  both  lasing  transitions.  The  figure 
of  merit  becomes  the  sum  of  this  product  from  both  transitions. 
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III.  Material  Characterization 


Choosing  a  laser  crystal  for  the  two  wavelength  laser  radar  application  which 
has  the  optimum  combination  of  large  absorption  bandwidth,  large  stimulated  emis¬ 
sion  cross  section,  and  long  upper  state  lifetime  is  difficult,  in  part,  because  of  the 
large  number  of  possible  materials.  These  materials  are  the  result  of  an  ongoing  quest 
for  the  ultimate  host  material  for  neodymium  which  will  yield  a  combination  of  high 
pumping  efficiency,  high  gain,  thermal  stability,  and  good  optical  quality  (12:387). 
A  numerical  analysis  to  choose  the  best  crystal  from  among  those  available  must 
rely  on  the  published  data  for  the  material  parameters  for  each  crystal.  Unfortu¬ 
nately,  these  reported  values  either  vary  tremendously,  as  shown  in  Tables  3.1  -  3.5, 
or  the  values  have  not  been  published.  So,  the  only  way  to  draw  meaningful  conclu¬ 
sions  about  the  relative  performance  of  these  available  crystals  is  to  measure  their 
performance  in  an  identical  laser  setup. 

The  cross  sections  given  in  Tables  3.1  -  3.5  are  effective  laser  cross  sections. 
The  effective  stimulated-emission  cross  section  is  the  spectroscopic  cross  section  mul¬ 
tiplied  by  the  occupation  of  the  upper  laser  level  relative  to  the  entire  4F3/2  manifold 
(6:50). 

Table  3.1  shows  the  values  for  Nd:YAG.  This  crystal  is  probably  the  most  com¬ 
monly  used  laser  material  for  lasers  at  1.06/xm,  and  it  is  often  used  as  the  standard 
against  which  new  materials  are  compared.  Nd:YLF  is  attractive  as  a  laser  material 
because  it  has  a  long  upper  state  lifetime  (see  Table  3.2).  Nd:YV04  has  been  stud¬ 
ied  recently  because  of  its  large  peak  absorption  cross  section  and  broad  absorption 
bandwidth  in  the  pump  region,  as  well  as  its  large  stimulated  emission  cross  section 
at  both  1.06/xm  and  1.32/xm  (14:5546;  15:298).  The  drawbacks  of  Nd:YVC>4  are  the 
nonavailability  of  high  quality  crystals  and  the  short  upper  state  lifetime,  which  in¬ 
dicates  that  it  does  not  store  the  pump  energy  as  well  as  the  other  materials  with 
longer  lifetimes.  Nd:GdV04  (neodymium  gadolinium  orthovanadate)  is  a  relatively 
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Table  3.1  Published  spectroscopic  data  for  Nd:YAG  (neodymium  yttrium  alu¬ 
minum  garnet).  The  stimulated  emission  cross  sections  are  the  effective 
cross  sections,  which  take  into  account  the  Boltzmann  distribution  of  the 
electrons  within  the  two  Stark  sub-levels  of  the  4F 3/2  manifold. 


Reference 

%  doping 

Tf  fas) 

A  =  1.064/im 
gse  10“ 19  cm2 

A  =  1.34/im 
gse  10-19  cm2 

(  18:561) 

244 

3.4 

(29:3067) 

1.0 

245 

3.2 

(  8:572) 

230 

6.5 

(  11:1443) 

1 

244 

3.4 

0.9 

(  17:376) 

1.1 

2.1 

(  5:90) 

1 

220 

2.8 

(  14:5547) 

1 

2.8 

(  28:6375) 

1.1 

240 

2.8 

(30:1869) 

1.0 

272±  12 

(  31:1314) 

1.0 

237  ±10 

this  work 

1.0 

237 

Table  3.2  Published  spectroscopic  data  for  Nd:YLF  (neodymium  yttrium  lithium 
fluoride). 


Reference 
(  18:561) 
(  8:572) 

(  11:1443) 
(  28:6375) 
this  work 


%  doping 


1.1 


Tf  (MS) 
480 
480 
480 
500 
504 


A  =  1.047/zm 
i tse  10~19  cm2 
3/7 
3.2 
3.7 
1.9 


A  =  1.313  /on 
Gse  10~19  cm2 


0.6 


1.0 


Table  3.3  Published  spectroscopic  data  for  Nd:YV04  (neodymium  yttrium 
orthovanadate). 


Reference 

%  doping 

Tf  (/Xs) 

A  =  1.064/xm 
use  10~19  cm2 

A  =  1.3/xm 
ctse  10“ 19  cm2 

(  18:561) 

1.0 

98 

6.2 

(29:3067) 

1.0 

97 

9.84±  1.4 

(  15:298) 

7±2 

(  8:572) 

98 

20.1 

(  5:90) 

1 

100 

15.6 

(  28:6375) 

1.1 

100 

15.6 

(  13:43) 

115 

10.7 

Table  3.4  Published  spectroscopic  data  for  Nd:GdV04  (neodymium  gadolinium 
orthovanadate). 


Reference 

%  doping 

Tf  (/Xs) 

A  —  1.064/xm 
i tse  10-19  cm2 

A  =  1.3/xm 
&se  10-19  cm2 

(  17:376) 

1.2 

90 

7.6 

1.8 

(  5:90) 

1 

94 

7.6 

(  13:43) 

«100 

>7.6 

Table  3.5  Published  spectroscopic  data  for  Nd:SVAP  (neodymium  strontium 
fluorovanadate). 


Reference 

%  doping 

rf  (ns) 

A  =  1.064/xm 
use  10-19  cm2 

A  =  1.3/im 
&se  10-19  cm2 

(  5:90) 

0.4 

220 

5.0 

(  13:43) 

215 

5 

(  28:6375) 

0.38 

230 

5 

this  work 

216 

Table  3.6  Published  spectroscopic  data  for  Nd:SFAP  (neodymium  strontium  fluo- 
rapatite).  The  stimulated  emission  cross  sections  are  the  effective  cross 
sections. 


Tf  (ns) 

A  =  1.064/xm 
asE  10-19  cm2 

A  =  1.3/xm 
Ose  10-19  cm2 

298 

5.4 

2.4 

Reference  %  doping 
(  32:43)  ~ 


Table  3.7  Refractive  Indices  for  several  laser  crystal  materials  @  1064nm. 


n0 

ne 

Nd:YAG 

1.81633 

N.A. 

(35:318) 

Nd:YLF 

1.4481 

1.4704 

(6:60) 

Nd:SVAP 

1.824 

1.809 

(13:43) 

Nd:YV04 

1.958 

2.168 

(13:43) 

Nd:SFAP 

1.629 

1.624 

(13:43) 

new  material  which  is  receiving  a  lot  of  attention  recently  because  it  exhibits  sim¬ 
ilar  spectroscopic  properties  to  that  of  Nd:YV04  (16:1071;  10:999;  17:373).  What 
makes  it  even  more  interesting  is  that  high  quality  crystals  are  easier  to  manufacture 
because  of  the  closeness  in  size  between  the  neodymium  atom  and  the  gadolinium 
atom  it  displaces  in  the  crystal  lattice  (17:373). 

Nd:SVAP  and  Nd:SFAP  both  have  cross  sections  at  1.06/im  which  are  ap¬ 
proximately  twice  the  values  published  for  Nd:YAG.  The  lifetime  of  Nd:SVAP  is 
only  around  10%  below  the  values  published  for  Nd:YAG  whereas  the  lifetime  for 
Nd:SFAP  is  approximately  20%  higher.  The  values  for  Nd:SFAP  indicate  that  it  will 
have  much  higher  gain  than  Nd:YAG  for  a  similar  amount  of  pump  power. 

The  final  section  in  Chapter  II  outlined  a  general  figure  of  merit  for  choosing 
the  best  material  for  the  dual  wavelength  laser  application.  Table  3.8  on  page  3-5 
shows  the  calculation  of  this  metric  for  some  of  the  most  promising  laser  crystals. 
The  values  used  to  calculate  the  metric  for  each  material  were  drawn  from  the  most 
recently  published  values  in  Tables  3.1  -  3.5.  There  is  no  entry  in  Table  3.8  due  to 
the  lack  of  a  published  cross  section  for  the  4F3/2  -  4Ii3/2  transition.  Based  on  the 
information  in  Table  3.8,  Nd:SFAP  is  the  optimum  material  for  a  dual  wavelength 
laser  system.  This  conclusion  should  be  supported  with  appropriate  experimental 
measurements.  Also,  this  conclusion  disregards  issues  such  as  beam  quality,  thermal 
lensing,  and  other  material  and  system  related  issues. 


Table  3.8  Sum  of  the  asE^f  product  from  both  the  4F3/2  -  4In/2  and  4F3/2  -  4Ii3/2 
transitions. 


O SE(1.06fim)Tf  +  OsE(1.3fim)Tf 

(10-17/is  cm2) 

Nd:YAG 

10.5 

Nd:YLF 

20.6 

Nd:GdV04 

8.5 

Nd:YV04 

17.0 

Nd:SFAP 

23.2 

8.1  Absorption  measurements 

The  measurements  of  the  polarized  absorption  for  Nd:YAG,  Nd:YLF,  and 
Nd:SVAP  in  the  wavelength  region  of  interest  for  laser  diode  pumping  are  shown 
in  Figure  3.1.  Nd:YAG  is  an  optically  isotropic  crystal  and  so  there  is  no  polariza¬ 
tion  dependence  on  the  absorption.  Nd:YLF  and  Nd:SVAP  are  both  uniaxial  crystals 
and  exhibit  different  absorption  characteristics  for  incident  radiation  polarized  along 
the  optic  axis  (E||c)  versus  perpendicular  to  it  (EJ_c).  The  peak  absorption  wave¬ 
length  for  Nd:SVAP  occurs  at  808. 5nm  for  both  polarizations;  however,  the  peak 
absorption  is  stronger  for  E|]c.  The  peak  absorption  for  Nd:YLF  occurs  at  792nm 
with  the  electric  field  polarized  parallel  to  the  optic  axis,  however,  the  absorption 
at  792nm  for  EJ_c  is  relatively  small.  In  order  for  a  pump  source  to  be  efficient  at 
this  wavelength,  it  must  linearly  polarized,  and  the  polarization  must  be  carefully 
aligned  to  the  optic  axis  of  the  crystal. 

The  absorptance  shown  in  each  plot  in  Figure  3.1  is  the  total  fraction  of  incident 
pump  power  absorbed  as  a  function  of  wavelength.  Both  the  Nd:YAG  and  Nd:YLF 
sample  were  0.5cm  long  whereas  the  Nd:SVAP  sample  was  only  0.3cm  long.  Even 
with  this  shorter  length,  the  peak  absorptance  in  Nd:SVAP  is  comparable  to  the  peak 
in  Nd:YAG  and  Nd:YLF  indicating  that  it  has  a  very  high  absorption  cross  section 
at  the  pump  wavelength.  The  wavelength  accuracy  of  the  Perkin  Elmer  Lambda  9 
Spectrophotometer  which  measured  the  absorption  data  is  ±  0.8nm. 
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Absorptance  Absorptance  Absorptance 


Figure  3.1  Polarized  absorption  spectrum  for  Nd:YAG,Nd:YLF,  and  Nd:SVAP. 
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Keithley  Model  427  LeCroy  9200A  precision 

Current  Amplifier  digital  oscilliscope 


Figure  3.2  Fluorescence  measurement  setup.  The  Ti:Sappire  pumped  laser  crys¬ 
tal  was  replaced  with  a  black  body  source  for  the  measurement  of  the 
relative  wavelength  response. 


3.2  Fluorescence  measurements 

Using  the  peak  absorption  wavelength  for  each  laser  crystal  as  indicated  in 
Figure  3.1  for  a  pump  source,  the  polarization-dependent  fluorescence  was  analyzed. 
The  measurements  were  made  with  the  experimental  setup  shown  in  Figure  3.2  on 
page  3-7.  The  0.5  meter  Jarrel-Ash  monochromator  contained  a  grating  blazed  for 
1.37/im.  The  measurement  results  are  shown  in  Figure  3.3  and  3.4  and  they  indicate 
the  wavelengths  where  each  crystal  will  lase.  These  results  also  give  an  indication 
about  the  number  of  transitions  that  are  possible  and  their  relative  strength. 

The  multiple  wavelengths  present  in  the  Nd:YAG  and  Nd:YLF  emission  are  in 
part  due  to  the  splitting  of  the  4F3/2  and  4In/2  multiplets.  This  splitting  is  present  in 
Nd:SVAP,  but  the  Stark  splitting  is  unusually  large  and  so  the  emission  is  dominated 
by  a  single  transition.  This  transition  is  stronger  than  in  Nd:YAG  because  it  occurs 
from  the  heavily  populated  lower  level.  This  lower  level  is  more  heavily  populated 
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than  the  equivalent  levels  in  Nd:YAG  due  to  the  Boltzmann  distribution  of  the 
excited  ions  in  the  4F3/2  multiplet  (13:44;  17:376). 

The  fluorescence  spectra  have  been  corrected  for  the  relative  wavelength  re¬ 
sponse  of  the  measurement  setup.  The  wavelength  response  was  measured  by  placing 
a  calibrated  blackbody  source  at  a  temperature  of  1000°C  in  front  of  the  aperture. 
This  source  was  directed  straight  into  the  monochromator  with  no  reflections.  The 
system’s  relative  response,  shown  in  Figure  3.5  on  page  3-11,  was  calculated  by 
dividing  the  measurement  of  the  output  of  the  blackbody  source  by  the  theoretical 
radiance  from  a  blackbody.  The  spectrum  for  each  material  was  essentially  corrected 
by  dividing  the  measured  spectrum  by  the  relative  response. 

Since  there  was  no  detector  available  which  could  measure  both  the  1.32/xm 
emission  from  the  laser  crystal  and  the  output  from  a  calibration  lamp,  calibration 
of  the  wavelength  scale  in  the  emission  measurements  was  a  two-step  process.  The 
monochromator  counter  reading  was  calibrated  using  a  neon  lamp.  The  location 
of  the  lines  in  the  second  and  third  order  of  the  grating  were  used  to  calibrate  the 
counter  in  the  vicinity  of  1.06/wn  and  1.32/xm.  During  a  spectrum  measurement, 
the  output  of  the  detector  was  amplified  using  a  current  amplifier  and  then  recorded 
using  the  LeCroy  9200A  oscilloscope.  This  recorded  output  was  matched  to  the 
monochromator  dial  setting  by  passing  a  bright  light  over  the  detector  when  the 
monochromator  dial  passed  through  specific  points  before  and  after  the  wavelength 
regions  of  interest.  The  measured  detector  trace  then  contained  two  spikes  which 
corresponded  to  specific  dial  settings  on  the  monochromator  allowing  the  calculation 
of  dial  settings  for  each  point  in  the  trace.  The  dial  settings  were  then  converted  to 
actual  wavelengths  from  the  earlier  calibration  using  a  neon  lamp. 

3.3  Lifetime  measurements 

The  fluorescence  lifetime,  the  lifetime  of  the  4F3/2  energy  state,  was  measured 
in  Nd:YAG,  Nd:YLF,  and  Nd:SVAP  with  the  setup  shown  in  Figure  3.6.  The  results 
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ed  units) 


Figure  3.3  Fluorescence  spectra  for  Nd:YAG,  Nd:YLF,  and  Nd:SVAP  in  the  vicin 
ity  of  1.06/zm. 
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Figure  3.4  Fluorescence  spectra  for  Nd:YAG,  Nd:YLF,  and  Nd:SVAP  in  the  vicin¬ 
ity  of  1.32/mi. 


u  1  1.05  1.1  1.15  1.2  1.25  1.3  1.35 
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Figure  3.5  Relative  system  response  of  the  setup  used  to  measure  the  fluorescence 
spectra  from  each  laser  crystal. 

are  plotted  on  a  log  scale  in  Figures  3. 7-3. 9.  The  measured  lifetimes  for  Nd:YAG 
and  Nd:SVAP  are  within  1.5%  of  the  published  values.  The  measured  lifetime  for 
Nd:YLF  which  was  504/zs,  however,  is  5%  higher  than  the  typical  published  value 
of  480/is  but  is  still  under  the  value  of  520/is  quoted  by  Dallas  as  typical  (28:6374) 
(I  have  not  found  any  other  publications  which  support  his  value).  One  possible 
explanation  for  this  larger  measured  value  is  radiation  trapping.  The  lifetime  of 
4Iii/2  energy  level  in  Nd:YLF  is  20ns,  three  orders  of  magnitude  larger  than  the 
225ps  in  Nd:YAG.  The  population  residing  in  this  lower  energy  state  is  available  for 
reabsorption  -  re-emission  causing  the  lifetime  to  appear  longer  than  it  actually  is. 

The  measured  decay  signal  was  corrected  by  subtracting  the  zero-offset  bias 
from  the  detector  signal  and  finding  the  best  straight  line  to  the  log  of  the  intensity 
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Figure  3.6  Experimental  setup  for  measuring  the  fluorescence  lifetime  of  Nd:YAG, 
Nd:YLF,  and  Nd:SVAP 

versus  time.  With  the  pump  beam  blocked,  the  intensity  is  assumed  to  decay  as: 

I(t)  =  /0exp  (3>1) 

where  Tj  is  the  fluorescence  lifetime,  the  lifetime  of  the  upper  population  level. 
Taking  the  natural  logarithm  of  the  above  equation  yields  a  straight  line  with  a 
slope  equal  to  the  fluorescence  lifetime.  The  displayed  curve  for  Nd:YAG  is  the 
result  of  averaging  100  traces  on  the  digital  oscilloscope;  the  curves  for  Nd:YLF  and 
Nd:SVAP  are  the  average  of  200  traces.  The  error  from  the  curve  fit  was  less  than 
1.  /zs  for  each  lifetime  measurement. 
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Intensity  (arb.  units) 


1 


Figure  3.7  Lifetime  measurement  of  Nd:YAG.  Both  best  fit  as  well  as  the  measured 
decay  curve  are  shown  (they  overlap).  The  displayed  data  is  the  average 
of  100  traces  taken  with  the  oscilloscope.  The  chopper  speed  used  in 
this  measurement  was  2394  Hz. 
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Intensity  (arb.  units) 


Intensity  (arb.  units) 


Figure  3.9  Lifetime  measurement  of  Nd:SVAP.  Both  best  fit  as  well  as  the  measured 
decay  curve  are  shown  (they  overlap).  The  displayed  data  is  the  average 
of  200  traces  taken  with  the  oscilloscope.  The  chopper  speed  used  in 
this  measurement  was  2112  Hz. 
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IV.  Experimental  results  &  analysis 

4-1  Experimental  setup 

A  perspective  view  of  the  optical  setup  is  shown  in  Figure  4.1  and  it  shows 
the  Burleigh  Wavemeter  used  to  monitor  the  wavelength  of  the  Ti:Sapphire  pump 
beam.  This  setup  proved  extremely  convenient  for  accurately  tuning  the  wavelength 
of  the  Ti:Sapphire  to  the  wavelength  of  peak  absorption  in  each  crystal.  The  neutral 
density  filter  wheel  was  used  to  vary  the  incident  pump  power  into  the  cavity.  The 
alignment  of  the  Ti'.Sapphire  pump  beam  and  the  calculation  of  the  position  of  the 
30cm  focusing  lens  is  covered  in  Appendix  A.  The  dichroic  beamsplitter  in  the 
“Y”  cavity  configuration  is  used  to  separate  the  two  different  wavelengths.  Since 
the  dichroic  beamsplitter  reflects  the  1.3/xm,  the  1.3/im  cavity  is  “L”  shaped.  The 
dichroic  beamsplitter  allows  the  1.06/im  beam  to  travel  through  so  the  1.06/im  cavity 
is  straight. 

Using  thermo-electric  coolers,  the  temperature  of  the  laser  crystal  was  main¬ 
tained  at  19.5±0.5°C  for  all  of  the  laser  materials.  Due  to  the  high  humidity  in  the 
laboratory,  water  started  to  condense  on  the  mount  if  the  temperature  was  dropped 
more  than  5°C  below  room  temperature.  This  prohibited  any  investigations  into  the 
temperature  sensitivity  of  the  laser  performance. 

4-1.1  Nd:YAG  setup.  The  Nd:YAG  crystal  was  a  cylindrical  rod  4mm  in 
diameter  and  5mm  long.  The  absorption  curve  in  Figure  3.1  indicated  that  the  best 
pump  wavelength  was  808nm.  The  Ti:Sapphire  was  tuned  to  this  wavelength  and  fo¬ 
cused  down  into  the  laser  rod.  After  the  cavity  was  lasing  at  1.06/im,  the  wavelength 
was  tuned  to  the  point  where  maximum  power  was  achieved.  This  wavelength  was 
found  to  be  808.54±0.6nm  as  measured  with  the  Burleigh  wavemeter.  When  setting 
up  for  Q-switched  operation,  the  optimum  angle  of  the  dichroic  surface  normal  with 
the  incident  beam  was  found  to  be  41.8±0.3°  for  maximum  1.06/im  output  power. 
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Figure  4.1  Optical  setup  for  the  laser  system. 

4.1.2  Nd:YLF  setup.  The  setup  for  Nd:YLF  proved  to  be  a  bit  more 
involved  than  for  Nd:YAG.  The  Nd:YLF  was  also  a  cylindrical  rod  4mm  in  diameter 
and  5mm  long.  Nd:YLF,  however,  is  a  uniaxial  crystal  and  this  particular  crystal 
was  an  “‘a-cut’”  crystal  meaning  that  both  polarizations  could  be  accessed.  The 
output  from  the  Ti:Sapphire  was  linearly  polarized  with  the  polarization  horizontal 
to  the  surface  of  the  optical  table.  The  first  step  in  setting  up  the  cavity  for  Nd:YLF 
was  to  tune  the  wavelength  of  the  Ti:Sapphire  to  the  peak  absorption  wavelength  of 
792nm,  as  indicated  in  Figure  3.1.  This  beam  was  focused  down  into  the  laser  rod 
using  the  30cm  focal  length  lens.  A  protected  silver  mirror  was  placed  on  the  other 
side  of  the  rod  to  reflect  the  residual  pump  beam  into  a  power  meter.  The  rod  was 
then  rotated  until  this  transmitted  power  was  minimized.  The  mirror  was  removed 
and  the  cavity  was  aligned.  Once  the  1.047/zm  cavity  was  lasing,  the  wavelength  of 
the  Ti'.Sapphire  was  adjusted  to  achieve  maximum  power  out  of  the  Nd:YLF  laser. 
When  setting  up  for  Q-switched  operation,  the  optimum  angle  of  the  dichroic  surface 
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normal  with  the  incoming  beam  for  Nd:YLF  was  found  to  be  45.0±0.5°  for  maximum 
1.047/im  output  power. 

4-1.3  Nd:SVAP  setup.  The  Nd:SVAP  crystal  was  a  3mm  cube  and  it  was 
also  an  ‘“a-cut”  ’  uniaxial  crystal.  The  brass  mount  for  this  crystal  was  made  with 
the  sample  5mm  too  low.  This  meant  that  the  top  of  the  high  reflector  needed 
to  be  tilted  towards  the  crystal,  the  result  was  that  the  pump  beam  had  a  higher 
transmission  through  the  high  reflector  for  this  location  and  tilt  angle. 

This  crystal  underwent  catastrophic  failure  during  the  setup  for  the  Q-switch 
measurements.  The  amount  of  incident  pump  power  was  approximately  1.0  Watt 
and  the  beam  diameter  in  the  crystal  was  270  ±  10  pm.  Visually,  the  damage 
appeared  to  be  faint  fracture  lines  internal  to  the  crystal. 

The  failure  of  the  Nd:SVAP  crystal  brings  to  light  one  aspect  of  the  laser  mate¬ 
rial  which  had  not  been  investigated  -  manufacturability  and  durability.  The  crystal 
itself  was  originally  to  be  a  4mm  by  4mm  crystal  but,  because  of  imperfections,  was 
cut  down  to  a  3mm  cube.  Manufacturers  of  this  type  of  crystal  have  found  that  they 
tend  to  fracture  during  the  cutting  process  so  they  tend  to  have  a  poor  yield  which 
contributes  to  a  relatively  high  cost(12:51). 

4-2  Optical  component  characterization 

This  next  section  steps  through  each  component  and  presents  the  important 
transmission  measurements  for  the  cavity  components  and  diagnostic  optics.  To  un¬ 
derstand  and  accurately  predict  the  behavior  of  a  laser  system,  the  characteristics  of 
the  components  which  make  up  that  system  must  be  known.  As  seen  in  Eq  (2.13), 
predicting  the  pump  power  necessary  to  achieve  lasing  threshold  requires  two  impor¬ 
tant  values;  the  amount  of  loss  that  each  component  causes  at  the  laser  wavelength, 
and  the  percentage  of  the  pump  power  incident  onto  the  cavity  which  is  actually 
absorbed  in  the  laser  gain  medium. 
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Table  4.1  Transmission  of  the  cavity  high  reflector  for  the  different  wavelengths. 


Laser  crystal 

Thr 

Wavelength 

Transmission 

(nm) 

Nd:YAG 

0.74  ±0.01 

808. lnm  0°  angle  of  incidence 

Nd:YLF 

0.791  ±  0.006 

792. lnm 

Nd:SVAP 

0.821  ±0.001 

808. 5nm  20°  angle  of  incidence 

4-2.1  High  reflector.  For  all  of  the  laser  measurements,  the  power  in  the 
incident  pump  beam  was  measured  after  the  30cm  focusing  lens.  Traveling  to  the 
right  toward  the  laser  crystal,  the  pump  beam  initially  passes  through  the  high 
reflector.  The  transmission  measurements  of  the  high  reflector  tjur  at  the  pump 
wavelengths  for  the  various  laser  crystals  are  shown  in  Table  4.1. 

The  high  reflector  is  also  the  first  element  in  both  laser  cavities.  The  reflectivity 
of  this  element  at  1.06/jm  was  measured  to  be 

THr  =  0.9998  ±  0.0001 

where  the  reflectivity  T  was  calculated  from  the  transmission  using  T  =  1  —  T 
(this  assumes  no  absorption  loss).  The  uncertainty  was  derived  from  the  standard 
deviation  of  multiple  measurements. 

4-2.2  Laser  Crystal.  After  passing  through  the  high  reflector,  the  pump 
beam  is  incident  onto  the  front  face  of  the  laser  crystal.  The  amount  of  pump 
power  which  is  reflected  from  this  surface  was  measured  using  the  setup  shown  in 
Figure  4.2.  The  uncertainties  in  the  reflectivity  values  came  from  the  standard 
deviation  of  multiple  measurements.  The  reflectivity  measurements  of  the  different 
crystals  are  shown  in  Table  4.2. 

After  passing  through  the  bulk  of  the  laser  crystal,  the  pump  beam  exits  the 
back  face,  which  is  assumed  to  have  the  same  power  reflectivity  as  the  front  face.  Two 
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Figure  4.2  Setup  for  the  transmission  measurements  of  the  laser  crystals. 


Table  4.2  Reflectivity  of  laser  crystals  at  the  appropriate  pump  wavelength. 


Laser  crystal 

Tic 

Reflectivity 

Wavelength 

(nm) 

Nd:YAG 

0.0498  ±  0.0002 

808. lnm 

Nd:YLF 

0.0050  ±  0.0008 

792. lnm 

Nd:SVAP 

0.042  ±  0.001 

808. 5nm 
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Figure  4.3  Setup  for  the  absorptance  measurements  of  the  laser  crystals. 


values  of  power  were  used  to  calculate  the  transmission  at  the  pump  wavelength.  The 
locations  where  they  were  measured  are  shown  in  Figure  4.3.  The  power  measured 
after  the  crystal  must  be  corrected  for  the  percentage  reflected  off  of  the  back  face  in 
order  to  accurately  calculate,  r)Abs.  The  fraction  of  incident  power  actually  absorbed 
in  the  crystal  is  given  by 


incident 


(1  —  Tic)  ~ 


VAbs  — 


*  transmitted 

(1  —  Tic ) 


Pi 


incident 


The  average  absorptance  of  each  laser  crystal  was  measured  at  its  pump  wave¬ 
length  using  the  experimental  setup  shown  in  Figure  4.3  on  page  4-6.  The  results  are 
shown  in  Table  4.3  on  page  4-7.  The  expression  for  calculating  the  amount  of  pump 
power  absorbed  in  the  laser  crystal  from  the  value  of  the  pump  power  measured 
before  the  high  reflector  becomes 


PAbsorbed  —  VAbs  VH  R  Pincident 


(4.1) 
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Table  4.3  Absorptance  of  laser  crystals  at  the  appropriate  pump  wavelength. 


Laser  crystal 

VAbs 

Absorptance 

Wavelength 

Nd:YAG 

0.88  ±  0.01 

808.  Inm 

Nd:YLF 

0.996  ±  0.001 

792.  lnm 

Nd:SVAP 

0.923  ±  0.009 

808. 5nm 

4-2.3  Dichroic  beamsplitter.  The  transmission  of  the  dichroic  beamsplitter 
measured  at  1.064/im  using  the  unpolarized  output  of  the  Nd:YAG  laser  was 

TDichroic(Nd:YAG)  =  0-90  ±  0.02 

The  amount  of  power  reflected  by  the  dichroic  was  measured  while  it  was  lasing 
in  the  Nd:YLF  cavity  at  the  optimum  angle.  The  amount  of  1.047/im  output  power 
was  back  propogated  into  the  cavity  to  calculate  the  amount  of  power  traveling  to  the 
left  and  right  in  the  cavity  just  on  the  other  side  of  the  output  coupler.  The  amount 
of  power  reflected  by  the  output  coupler  was  calculated  from  the  reflectivity  of  the 
output  coupler  and  the  output  power.  This  reflected  power  was  used  as  the  incident 
value  onto  the  dichroic  beamsplitter.  The  transmission,  calculated  by  comparing 
this  value  with  the  power  reflected  by  the  dichroic  beamsplitter,  was 


TDichroic(Nd:YLF)  =  0.982  ±  0.002 

The  wavelength  and  polarization  dependent  transmission  of  the  dichroic  beam¬ 
splitter  was  also  measured  using  a  Perkin  Elmer  Lambda  9  Spectrophotometer. 
These  results  are  shown  for  1.06/im  in  Figure  4.4.  These  curves  tended  to  shift 
horizontally  in  wavelength  as  the  incident  angle  was  changed.  This  effect  caused  the 
dichroic  beamsplitter  to  have  a  different  optimum  angle  in  the  cavity  than  Nd:YLF. 
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Figure  4.4  Transmission  of  the  dichroic  beamsplitter  in  the  vicinity  of  1.06/im  at 
45°  angle  of  incidence.  These  curves  demonstrate  the  relative  transmis¬ 
sion  levels  at  different  wavelengths  and  were  taken  with  a  Perkin  Elmer 
Lambda  9  UV/VIS/NIR  Spectrophotometer. 

The  transmission  was  also  measured  in  the  vicinity  of  1.3/zm  with  the  values  con¬ 
verted  to  a  reflectivity  using  F=l-T;  the  results  are  shown  in  Figure  4.5  on  page 
4-9. 


4-2-4  Q-switches.  The  transmission  of  the  1.06/xm  Q-switch  was  measured 
with  the  output  of  the  Nd:YAG  laser  and  found  to  be 


TW6QS  =  0.95  ±  0.02 
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Percent  Reflectivity 


Figure  4.5  Reflectivity  of  the  dichroic  beamsplitter  in  the  vicinity  of  1.32/im  at  45° 
angle  of  incidence.  These  measurements  were  made  with  a  Perkin  Elmer 
Lambda  9  UV/VIS/NIR  Spectrophotometer. 


Table  4.4  Transmission  of  the  long  pass  filter  at  the  laser  wavelengths 


A 

Transmission 

Tlp 

1.047/mi 

0.58  ±  0.01 

1.064/xm 

0.636  ±  0.006 

1.32/im 

0.856±  0.005 

The  transmission  of  both  Q-switches  was  also  measured  with  the  Perkin  Elmer 
Lambda  9  Spectrophotometer  and  the  results  were 


^io6QS  —  0-98  ±  0.5 


and 

^132 qs  =  0.98  ±  0.8 

4-2.5  Output  couplers.  The  reflectivities  of  the  two  1.06/xm  output  couplers 
were  measured  as 


r95%  =  0.95  ±  0.01 


r70%  =  0.700  ±  0.008 


4-2.6  Diagnostic  optics.  To  filter  out  the  pump  wavelength  the  ouput  of 
each  laser  was  passed  through  a  long  pass  filter  made  from  Schott  RG-1000  glass. 
The  transmission  of  this  optic  at  the  laser  wavelengths  was  measured  and  the  results 
are  shown  in  Table  4.4. 
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The  reflectivity  of  the  protected  silver  fold  mirrors,  used  to  direct  the  laser 
beam  into  the  power  meter,  was  measured  at  an  incident  angle  of  45°  to  be 

Tfoid  =  0.953  ±0.005 

where  the  uncertainty  is  derived  from  multiple  measurements. 

4.3  CW  laser  Performance 

This  section  presents  the  results  from  measurements  made  with  the  laser  cavity 
lasing  in  CW  mode.  The  results  are  used  to  verify  the  amount  of  loss  caused  by  each 
element  in  the  laser  cavity  and  to  estimate  the  pumping  efficiency.  Both  the  pumping 
efficiency  and  the  cavity  loss  values  are  necessary  to  accurately  model  the  Q-switched 
laser  performance. 

The  first  set  of  CW  measurements  were  made  with  only  the  laser  crystal  in  the 
laser  cavity.  These  measurements  can  be  used  to  estimate  the  combined  loss  from 
the  high  reflector  and  from  the  laser  crystal.  Eq  (2.14)  on  page  2-7  gives  an  equation 
for  determining  the  cavity  loss  from  measurements  of  the  laser  pump  threshold  for 
different  output  coupler  reflectivities.  Two  output  coupler  reflectivities,  95%  and 
70%  were  available  for  1.06/xm  measurements  but  only  one,  96%,  at  1.3/rni.  Thresh¬ 
old  measurement  results  using  the  available  output  couplers  are  shown  in  Table  4.5. 
The  wavelength  measurements  at  1.06/im  were  made  with  the  Burleigh  Waveme- 
ter,  and  the  wavlength  measurements  at  1.32/zm  were  made  with  the  Jarrel-Ash 
monochromator.  The  1.3/un  CW  output  of  both  Nd:YAG  and  Nd:YLF  contained 
two  wavelengths,  which  is  common  for  Nd-based  lasers  operating  in  the  1.3/mi  re¬ 
gion  (34:52).  Figure  3.4  on  page  3-10  shows  the  measured  emission  spectra  for  these 
materials.  Both  materials  have  two  dominant  energy  level  transitions  that  are  on 
the  same  order  of  magnitude. 
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Table  4.5  CW  Laser  threshold  and  wavelength  data 


Output  Coupler 
(%  Reflectivity) 

Xp 

(nm) 

Material 

(/zm) 

P  Threshold 

(mW) 

95 

808.1 

Nd:YAG 

1.064 

109±3 

70 

808.1 

Nd:YAG 

1.064 

670±20 

96 

808.1 

Nd:YAG 

1.338/1.318 

490±10 

95 

792.1 

Nd:YLF 

1.047 

78±2 

70 

792.1 

Nd:YLF 

1.047 

620±20 

96 

792.1 

Nd:YLF 

1.313/1.321 

270±10 

95 

808.5 

Nd:SVAP 

1.065 

82±2 

96 

808.5 

Nd:SVAP 

1.337 

280±10 

Table  4.6  Estimates  for  2r\e^iCitncylgo  sE^f  and  residual  cavity  loss  (L) 


Material 

L 

Nd:YAG 

0.545 

0.011  ±  0.003 

Nd:YLF 

0.561 

-0.005  ±  0.002 

The  two  threshold  values  at  1.06/zm  for  Nd:YAG  and  at  1.047/xm  for  Nd:YLF 
were  used  in  Eq  (2.14)  to  determine  a  value  for  the  total  cavity  loss  for  each  crys¬ 
tal.  The  results  of  these  calculations  are  shown  in  Table  4.6.  The  total  cavity  loss 
(excluding  the  output  coupler)  for  both  crystals  is  extremely  small.  The  value  for 
Nd:YLF  is  less  than  zero,  which  implies  gain.  Since  this  is  not  possible,  there  must 
be  some  uncertainty  in  the  threshold  pump  powers  that  is  not  evident  from  the 
measurements. 

The  next  set  of  CW  measurements  were  power  conversion  measurements  (out¬ 
put  laser  power  versus  input  pump  power),  and  these  were  also  made  using  the  simple 
cavity  configuration.  The  data  from  these  measurements  is  plotted  in  Figures  4.6 
and  4.7.  As  expected,  the  increased  fluorescence  lifetime  in  Nd:YLF  produces  a  lower 
threshold  and  higher  output  for  a  given  pump  power  than  Nd:YAG  -  but  not  as  high 
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as  would  be  predicted  from  a  material  with  over  twice  the  fluorescence  lifetime.  The 
poor  slope  efficiency  performance  of  NdtSVAP  at  1.06//m  is  unexpected  because  it 
has  a  higher  stimulated  emission  cross  section-fluorescence  lifetime  product  for  the 
1.06/im  transition  than  Nd:YAG  but  less  than  that  of  Nd:YLF.  This  relative  per¬ 
formance  is  reflected  in  the  1.3/im  results.  Although  there  is  no  published  cross 
section  for  Nd:SVAP  at  1.3/wn,  based  on  the  results  shown  in  Figure  4.7  and  the 
threshold  measurement  results  in  Table  4.5,  the  cross  section  should  be  much  larger 
than  both  Nd:YAG  and  Nd:YLF.  Comparing  the  threshold  pump  power  of  Nd:SVAP 
with  Nd:YLF  from  Table  4.5,  shows  the  threshold  for  NdtSVAP  is  only  4mW  higher 
in  the  identical  cavity.  Since  the  threshold  pump  power  is  inversely  proportional  to 
the  product  of  the  stimulated  emission  cross  section  and  fluorescence  lifetime  and 
since  the  lifetime  of  NdtYLF  is  over  twice  that  of  NdtSVAP,  the  conclusion  is  that 
NdtSVAP  has  a  much  higher  cross  section. 

The  1.05/zm  output  of  the  NdtYLF  crystal  was  found  to  be  polarized  at  least 
1000:1,  which  was  based  on  measurements  with  a  calcite  polarizer.  This  polarization 
was  along  the  optic  axis  of  the  crystal  oriented  horizontal  to  the  table.  This  orien¬ 
tation  was  required  in  order  to  be  aligned  with  the  polarization  of  the  TitSapphire 
pump  laser  (as  described  in  subsection  4.1.2). 

The  last  set  of  CW  measurements  focus  on  only  the  NdtYAG  cavity  and  were 
designed  to  provide  an  additional  measurement  of  the  intracavity  loss  caused  by  each 
component.  Table  4.7  shows  the  measured  slope  efficiencies  and  power  thresholds 
when  different  components  are  placed  in  the  cavity.  These  values  correspond  to  the 
laser  output  power  vs  absorbed  pump  power  curves  shown  in  Figure  4.8  on  page 
4-17.  The  pump  threshold  values  were  calculated  from  the  curves  by  extrapolating 
them  to  the  horizontal  axis.  Figure  4.8  shows  the  increase  in  threshold  pump  power 
as  each  element  is  added  to  the  cavity.  The  cavity  containing  both  the  Q-switch 
and  the  dichroic  beamsplitter  is  identical  to  the  one  used  to  make  the  Q-switch 
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Output  Power  (Watts) 


Figure  4.6  Power  conversion  measurements  for  1.06/mi.  The  same  cavity  was  used 
for  each  material. 
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Output  Power  (Watts) 


Absorbed  Pump  Power  (Watts) 

Figure  4.7  Power  conversion  measurements  for  1.32/xm.  The  same  cavity  was  used 
for  each  material. 
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measurements,  and  so  a  measurement  of  the  loss  in  this  cavity  is  useful  in  predicting 
the  Q-switched  laser  performance. 

The  equation  for  the  CW  laser  threshold  in  Eq  (2.13)  can  be  used  to  determine 
the  unknown  loss  caused  by  a  cavity  element.  The  equation  for  the  threshold  pump 
power  when  all  of  the  cavity  losses  are  known  is  given  by 


PThreshi 


-V0ohl>Lhl  (^J^2^lT^nowr^ 


G  S  E^fVef  ficiency  l  g 


(4.2) 


where  Tkuowu  is  the  known  one  way  transmission  of  the  components  in  the  cavity. 
The  threshold  pump  power  for  a  cavity  containing  an  unknown  loss  is 


PThreshi  — 


— Vqq hVL  In  (  ^/r2  r  1  PjjnWwn^Knowri) 

& S  E^fVef  ficiencylg 


(4.3) 


where  Tunknown  is  the  unknown  one  way  transmission  of  the  added  component  (s). 
Taking  the  ratio  of  these  two  equations  and  solving  for  Tunknown  yields 


exp 

Tunknown  = 


Threshi 

o 

Threshi 


in  (v^rvnuT)  } 


v/riiV72 


Known 


(4.4) 


This  expression  was  used  to  calculate  the  transmission  values  of  the  different  com¬ 
ponents  presented  in  Table  4.7.  The  value  for  the  transmission  of  the  dichroic  cal¬ 
culated  from  the  lasing  threshold  of  the  Nd:YAG  laser  is  within  the  experimental 
uncertainty  of  the  value  measured  with  the  unpolarized  Nd:YAG  output  as  a  source 
(see  subsection  4.2.3  on  page  4-6).  The  value  for  the  transmission  of  the  1.06/zm 
Q-switch  calculated  from  the  threshold  change  (Table  4.7)  is  4%  higher  than  the 
value  measured  in  subsection  4.2.4.  There  are  two  possible  explanations.  The  first 
is  a  problem  that  plagued  measurements  throughout  this  work,  and  that  is  contam¬ 
ination.  When  the  beam  sizes  are  on  the  order  of  less  than  a  millimeter,  the  dust 
particles  in  the  air  have  a  significant  impact  on  power  measurements.  A  small  dust 
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Figure  4.8  Power  conversion  measurements  for  NdrYAG  at  1.06/zm  with  different 
combinations  of  components  in  the  cavity. 

particle  on  either  face  of  the  dichroic  could  significantly  effect  the  transmission  and 
cause  it  to  be  lower  than  actual.  The  second  possible  explanation  is  that  by  placing 
the  Q-switch  into  the  cavity,  with  its  refractive  index  on  the  order  of  1.68,  it  caused 
the  effective  cavity  length  to  decrease.  Decreasing  the  effective  cavity  length  causes 
the  cavity  to  become  more  stable,  it  also  increases  the  effective  spot  size  of  the  laser 
mode  in  the  laser  crystal.  If  the  pump  beam  size  was  larger  than  the  laser  mode 
in  the  laser  crystal,  increasing  the  laser  mode  size  would  increase  the  beam  overlap 
efficiency  according  to  Eq  (2.10)  on  page  2-6. 

Figures  4.9  and  4.10  show  the  change  in  the  output  power  curve  for  Nd:YAG 
and  Nd:YLF,  respectively,  when  the  output  coupler  reflectivity  is  changed  from  70% 
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Table  4.7  Nd:YAG  laser  threshold  measurements  with  different  combinations  of 
components  in  the  cavity. 


Components 
in  the  cavity 

Wavelength 

(/mi) 

Slope 

Efficiency 

% 

Threshold 

Pump  Power 
(W) 

Transmission 
of  element  (s) 

% 

Laser  crystal  only 

1.06 

36. 4±  0.09 

0.1041  ±  0.0002 

99.2  (Table  4.6) 

Dichroic 

1.06 

29.0  ±1 

0.40  ±  0.006 

91.9 

Q-switch 

1.06 

36.4  ±  0.8 

0.134  ±  0.001 

99.2 

Q-switch  and  Dichroic 

1.06 

33  ±  4 

0.502  ±  0.008 

89.3 

Dichroic 

1.32 

81  ±  3. 

0.640  ±  0.001 

Q-switch 

1.32 

75  ±  0.2 

0.699  ±  0.002 

99.8 

to  95%.  Examining  the  curves,  the  slope  efficiency  does  not  change  between  the 
output  coupler  values.  The  reason  for  this  is  found  by  examining  the  equation  for 
the  slope  efficiency  (Eq  (2.9)),  repeated  here  for  convenience 


aS  —  T 


(1  r inefficiency 


1  +  yj ri TLmP  t1  \l^^2TComp\ 


If  the  greatest  loss  is  due  to  the  output  coupler  and  the  other  components  in  the 
cavity  have  very  little  loss,  that  is  their  collective  transmission  is  close  to  1,  then  the 
slope  efficiency  becomes 


_  (l  r2)  TJefficiency 

(l  +  VT^)  (l  —  yf^t) 

Combining  the  two  expressions  in  the  denominator,  this  becomes, 


as 


Vef  ficiency 
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Figure  4.9  Power  conversion  measurements  measurements  for  Nd:YAG  at  1.06//m 
for  two  different  output  couplers.  The  two  overlapping  curves  on  top  for 
the  95%  reflective  coupler  are  simply  repeated  measurements  and  they 
demonstrate  the  repeatability.  The  curve  to  the  right  was  measured 
using  a  70%  reflective  output  coupler. 

which  is  independent  of  the  output  coupler  reflectivity.  From  the  similar  slope  effi¬ 
ciencies  for  the  two  output  coupler  values  shown  in  Figures  4.9  and  4.10,  it  can  be 
concluded  that  the  cavity  used  to  make  these  measurements  was  very  low  loss. 


4-4  Single  wavelength  Q- switched  results 

This  section  outlines  a  set  of  Q-switched  measurements  made  for  Nd:YAG  and 
Nd:YLF  while  lasing  only  at  one  wavelength.  In  the  next  section,  the  model  output 
is  fit  to  the  measured  results  at  1.06/xm  for  Nd:YAG  using  the  results  from  the  CW 
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Output  Power  at  1 .047  microns  (Watts) 


Figure  4.10  Power  conversion  measurements  for  Nd:YLF  at  1.047/im  with  two  dif¬ 
ferent  output  coupler  reflectivities. 
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Figure  4.11  Measurement  setup  for  the  Q-switched  results. 

measurements  outlined  in  the  last  section,  the  energy  measurements  presented  in 
this  section,  and  the  theory  developed  in  Chapter  II.  The  lasing  results  for  Nd:YAG 
at  1.3 pan  and  for  Nd:YLF  at  both  1.047/xm  and  1.3  pan  are  merely  presented,  no 
attempt  was  made  to  fit  the  model  to  the  measured  results. 

Figure  4.11  shows  a  diagram  of  the  general  measurement  setup  used  for  all  of 
the  Q-switched  results.  The  1.06/xm  results  were  obtained  by  placing  a  beam  block 
in  the  1.3/xm  cavity  and  vice  versa.  The  measurements  for  Nd:YAG  at  1.064/mi  are 


4-21 


Table  4.8  Q-switched  Nd:YAG  single  wavelength  performance  at  1.06/rm. 


Incident 

Pump 

Power 

(W) 

Absorbed 

Pump 

Power 

(W) 

Pulse 

Repetition 

Rate 

(Hz) 

Pulse 

Energy 

Pulse 

Width 

(ns) 

1.45 

0.92±0.03 

500 

55  ±2 

36.4  ±0.2 

1.45 

0.92±0.03 

500 

56  ±2 

36.3±0.5 

1.30 

0.83±0.03 

500 

55  ±2 

40.1±0.8 

1.30 

0.83±0.03 

1000 

52  ±2 

42.1±0.6 

1.30 

0.83  ±0.03 

4000 

36  ±1 

57.1±0.8 

Table  4.9  Q-switched  Nd:YAG  single  wavelength  performance  at  1.32/rm. 


Incident 

Pump 

Power 

(W) 

Absorbed 

Pump 

Power 

(W) 

Pulse 

Repetition 

Rate 

(Hz) 

Pulse 

Energy 

OiJ) 

Pulse 

Width 

(ns) 

1.45 

0.92±0.03 

500 

45  ±2 

235  ±1 

1.45 

0.92±0.03 

500 

58  ±2 

228±4 

1.30 

0.83±0.03 

500 

37  ±1 

303±3 

1.30 

0.83±0.03 

1000 

35  ±1 

257±5 

1.30 

0.83  ±0.03 

4000 

21.8  ±0.6 

475±4 

shown  in  Table  4.8.  Figure  4.12  shows  a  plot  of  some  of  these  values  along  with  the 
results  of  the  model.  The  cavity  measurements  are  shown  in  Table  4.10. 

The  pulse  energies  were  calculated  from  an  average  power  measurement  of 
the  output.  There  is  no  commercially  available  energy  meter  which  can  measure 
the  amount  of  energy  in  each  pulse  at  repetition  rates  of  over  500Hz.  In  order  to 
calculate  an  average  pulse  energy,  the  average  power  of  the  laser  output  is  divided 
by  the  pulse  repetition  frequency. 
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Pulse  energy(microJoules)  Pulse 


Table  4.10  Cavity  configuration  for  the  Q-switched  Nd:YAG  data  -long  cavity. 


Pulse  Repetition  rate  (Hz) 


Figure  4.12  Q-switched  measurements  made  with  Nd:YAG  at  1.06  microns.  The 
model  results  are  plotted  with  the  solid  line  and  the  measured  results 
are  marked  with  and  “o” .  The  absorbed  pump  power  was  830mW. 
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Pulse  Energy 


‘wv500  1000  1500  2000  2500  3000  3500  4000 

Repetition  Rate  (Hz) 


Figure  4.13  Measurement  results  for  Nd:YAG  operating  only  at  1.34/im.  The  ab¬ 
sorbed  pump  power  was  830mW. 


The  Q-switched  measurements  of  the  Nd:YAG  1.32/xm  cavity  shown  in  Ta¬ 
ble  4.9  and  plotted  in  Figure  4.13.  Because  each  plot  has  only  three  points,  it  is 
difficult  to  draw  any  conclusions  with  confidence.  However,  the  pulse  energy  is  ex¬ 
pected  to  stay  constant  until  a  pulse  repetition  rate  equal  to  or  4.2kHz.  The  plots 
in  Figures  4.13  and  4.12  show  a  steady  decrease  in  energy  to  this  point  (no  data  was 
taken  beyond  this)  and  this  is  probably  due  to  the  effects  of  thermal  lensing  in  the 
laser  crystal. 

The  measured  results  from  Nd:YLF  shown  in  Tables  4.11  and  4.12  are  plotted 
in  Figures  4.14  4.15  for  the  1.05/im  and  1.32/im  cavities  respectively.  Similarly  to 
Nd:YAG,  the  pulse  energy  is  expected  to  roll  off  with  increasing  repetition  rate  at  a 
rate  corresponding  to  around  which  is  2.0kHz  for  Nd:YLF.  The  point  where  the 
roll  off  starts  to  occur  in  Figures  4.14  4.15  is  1.0kHz.  This  roll  off  does  match  the 
experimental  data  published  by  Grossman(3:623). 

Figures  4.14  and  4.15  also  show  a  decrease  in  pulse  energy  for  decreasing  rep¬ 
etition  rate  which  is  contrary  to  the  equations.  The  reason  for  this  decrease  can  be 
seen  by  examining  the  amount  of  energy  stored  in  the  material.  The  pump  beam 
was  focused  down  to  a  spot  size  of  135/im  inside  the  laser  crystal.  Previous  measure¬ 
ments  show  99.6%  of  the  incident  pump  power  absorbed  by  the  5mm  long  crystal. 
This  indicates  an  absorption  coefficient  0(Nd:YLF  of 


In  (1  —  .996)  , 

OiNd-.YLF  = - J -  =  11cm 

l9 


The  amount  of  power  absorbed  in  the  first  1  mm  of  the  crystal  is  given  by 


PAbsJlmm  —  Pine  (1  exp  (  (0.1cm)  Olpjd-.YLF )) 


where  P/nc  is  the  amount  of  power  in  the  incident  pump  beam.  This  can  be  converted 
to  the  energy  stored  in  the  4F3/2  state  per  unit  volume  using  Eqs  (2.2)  and  (2.3 


Wstored  =  p  (1  -  exp  (-o.lcma^y^))  Tf  (l  -  exp  (-^))  hv^ 

V  olimm  Vollmm,  hl/p 

where  t  is  the  population  build  up  time,  and  Volimm  is  the  volume  of  the  pumped 
region  in  the  first  1mm  of  the  crystal. 

Using  the  incident  pump  power  of  1.1 7W  used  to  generate  the  data  in  Fig¬ 
ures  4.14  and  4.15,  and  inserting  values  for  the  other  parameters  yields 


Wstored 

V  olimm 


11.5 


^1  —  exp 


Joules 

cm-3 


This  expression  is  a  function  of  the  pulse  repetition  rate.  The  effect  of  decreasing 
pulse  energy  begins  at  pulse  repetition  rate  of  1kHz  which  gives  a  stored  energy  of 


Wstored  _  g  g  Joules 
Volimm  '  cm"3 


This  significant  concentration  of  energy  causes  negative  thermal  lensing  in  Nd:YLF. 
This  negative  thermal  lensing  increases  the  effective  spot  size  of  the  laser  mode  in 
the  laser  crystal.  This  increase  in  the  spot  size  increases  the  threshold  population 
inversion  and  decreases  the  amount  of  extracted  energy  according  to  Eq  (2.22) 


4-5  Analysis  of  1.06yim  Nd:YAG  Q-switched  results 

The  model  used  to  analyze  these  results  was  the  same  one  used  to  design  the 
initial  system  and  is  shown  in  Appendix  B.  The  results  shown  in  Figure  4.12  on 
page  4-23  demonstrate  that  the  model  (solid  line)  gives  a  pulse  energy  and  a  pulse 
width  which  are  approximately  half  of  the  measured  values. 

Examining  the  energy  per  pulse  first,  a  minimum  value  for  the  amount  of 
absorbed  pump  power  that  is  needed  to  produce  these  energies  can  be  calculated 
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Table  4.11  Q-switched  Nd:YLF  single  wavelength  performance  at  1.047/im.  The 
data  for  the  cavity  configuration  is  given  in  Table  4.13 


Incident 

Pump 

Power 

(W) 

Absorbed 

Pump 

Power 

(W) 

Pulse 

Repetition 

Rate 

(Hz) 

Pulse 

Energy 

(pj) 

Pulse 

Width 

(ns) 

1.483 

1.17±0.04 

89.9 

87  ±4 

24.8  ±0.1 

1.483 

1.17±0.04 

150.6 

113±4 

24.3  ±0.2 

1.483 

1.17±0.04 

1061 

122  ±4 

26.1  ±0.6 

1.483 

1.17±0.04 

2080 

101  ±4 

28.4  ±0.3 

1.483 

1.17±0.04 

3000 

86  ±3 

32.6  ±0.6 

1.483 

1.17±0.04 

4000 

72  ±3 

35.9  ±0.3 

1.483 

1.17±0.04 

5000 

61  ±2 

39.8  ±0.3 

1.485 

1.17±0.04 

10680 

33  ±1 

60  ±1 

1.300 

1.02±0.03 

10680 

23.3  ±0.8 

78  ±2 

1.179 

0.93±0.03 

10680 

14.3  ±0.5 

123±4 

Table  4.12  Q-switched  Nd:YLF  single  wavelength  performance  at  1.32//m.  The 
data  for  the  cavity  configuration  is  given  in  Table  4.13. 


Incident 

Pump 

Power 

(W) 

Absorbed 

Pump 

Power 

(W) 

Pulse 

Repetition 

Rate 

(Hz) 

Pulse 

Energy 

w 

Pulse 

Width 

(ns) 

1.485 

1.17±0.04 

89.9 

47  ±2 

143  ±4 

1.485 

1.17±0.04 

150.6 

52±2 

134  ±3 

1.485 

1.17±0.04 

1009 

51  ±2 

156  ±5 

1.485 

1.17±0.04 

2000 

38  ±1 

204  ±4 

1.485 

1.17±0.04 

3000 

32  ±1 

233  ±7 

1.485 

1.17±0.04 

4000 

24.8  ±0.8 

265  ±8 

1.485 

1.17±0.04 

5000 

24.0  ±0.8 

300±10 

1.295 

1.02±0.03 

5000 

12.4  ±0.5 

422  ±9 
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Table  4.13  Cavity  configuration  for  the  Q-switched  Nd:YLF  data-long  cavity. 


li=1.0cm 

l4=6.7cm 

12=4. 2cm 

l5=1.9cm 

13=4. 1cm 

16=8. 5cm 

Output  coupler 

Cavity  length 

1.32//m 

T=4%,  20cmcc 

18.1cm 

1.06/zm 

T=30%,  20cmcc 

18.8cm 

Figure  4.14  Measured  results  for  Nd:YLF  operating  only  at  1.05/xm.  The  absorbed 
pump  power  was  1.17W. 
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Pulse  Energy  (microJoules) 


Figure  4.15  Measured  results  for  Nd:YLF  operating  only  at  1.32/mi.  The  absorbed 
pump  power  was  1.17W. 
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using  these  results  and  the  theory  developed  in  Chapter  II.  Table  4.8  on  page  4-22 
indicates  a  measured  pulse  energy  of  55±2/zJ  for  an  incident  pump  power  of  1.3W 
at  a  repetition  rate  of  500Hz.  The  expression  for  the  pulse  energy  given  in  Eq  (2.19) 
(repeated  here  for  convenience)  gives  the  pulse  energy  as  a  function  of  the  change  in 
population  inversion. 

Wout  =  ( rii  -nf)  hvL——— 

Qtotal 

Solving  this  equation  for  (nj  —  rif)  yields  the  change  in  the  population  inversion  that 
it  took  to  create  the  measured  pulse  energy: 


Q-totalW out 
Otexth'^L 


(4.5) 


Inserting  the  measured  energy  and  the  cavity  loss  from  Table  4.7  on  page  4-18  into 
Eq  (4.5)  gives  (nj  —  nj)  =  1.3  x  1014.  The  value  for  nth,  calculated  using  Eq  (2.7), 
is  nth  =  4.5  x  109.  Using  these  values  in  Eq  (2.22)  gives  the  fraction  of  population 
inversion  left  at  the  end  of  the  Q-switched  pulse 


rii 


which  says  that  n/  is  much  smaller  than  rii  and  so  it  is  safe  to  approximate  (nj  —  rif) 
as  just  rii  or 


m  &  ( rii  -  rif)  = 


Q 'total  ^ out 
01  ext  hl^L 


(4.6) 


This  value  can  be  used  in  the  expression  for  the  population  inversion  versus 
time  given  in  Eq  (2.2)  and  the  pump  rate  given  in  Eq  (2.3)  to  calculate  the  necessary 
absorbed  pump  power  to  produce  this  population  inversion.  This  absorbed  pump 
power  is 

P Pump  ’  Ve f  ficie.ncy  —  1.023  Watts 


This  is  greater  than  the  0.827  Watts  calculated  from  the  measured  transmission  val¬ 
ues.  This  difference  brings  to  light  one  difficulty  encountered  in  making  accurate 
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measurements  of  the  transmission  values  used  to  determine  the  amount  of  absorbed 
pump  power.  At  the  time  the  measurements  were  made,  the  only  power  meter  avail¬ 
able  with  the  necessary  dynamic  range  was  a  Coherent  Fieldmaster  with  a  model 
LM-2  silicon  detector  which  required  a  1000:1  attenuator  for  measurements  above 
50mW.  What  made  accurate,  repeatable  measurements  difficult  is  that  the  displayed 
power  varied  by  up  to  ±  30%  depending  on  the  position  and  angle  of  arrival  of  the 
beam  on  the  surface  of  the  attenuator.  This  places  an  uncertainty  on  the  measure¬ 
ments  which  were  made  absorptance  of  the  pump  beam  and  the  transmission  of  the 
optics  in  the  cavity,  and  leads  to  an  analysis  to  determine  their  relationship  from 
the  measurements  of  the  pulse  energy. 

Going  back  through  the  analysis,  the  two  most  uncertain  values  which  have  the 
greatest  impact  on  both  the  pulse  energy  and  the  pulse  length  are  the  total  internal 
cavity  losses,  atotai ,  and  the  absorption  efficiency  of  the  measured  pump  power, 
Veffidency  The  total  absorption  efficiency  includes  a  number  of  different  efficiency 
factors 

Veffidency  ~  TjjRVAbsVB 


Solving  Eq  (2.2),  Eq  (2.3),  and  Eq  (4.5)  for  these  two  values  yields  a  ratio  given  by 


V efficiency  out 

autcl  Vl  aeitT/  (l  -  <T  v  )  PPmp 


(4.7) 


_1 _ 

rep-rate 

for  single  wavelength  operation). 

Calculating  this  ratio  for  the  values  in  Table  4.8  yields 


where  t  is  the  amount  of  time  the  population  is  allowed  to  build  up  (^[[^r 


Veffidency  _  4  g  -j.  Q  2 
total 


Using  these  values  in  the  model  gives  the  results  shown  in  Figure  4.16.  The 
measured  pulse  energies  fall  within  the  bounds  of  the  ratio  but  the  calculated  pulse 
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Pulse  energy(microJoules)  Pulse  width  (ns) 


Figure  4.16  Model  results  for  pulse  width  and  pulse  energy  as  functions  of  pulse 
repetition  frequency.  The  multiple  curves  represent  the  variation  in  the 
value  of  the  ratio  of  absorption  efficiency  to  total  cavity  loss  calculated 
from  the  data. 
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width  is  still  approximately  half  of  the  measured  values.  One  possibility  is  that  the 
size  of  the  beam  in  the  laser  cavity  is  larger  than  the  TEMoo  spot  size  used  in  the 
model.  Using  the  nominal  value  for  the  ratio  and  increasing  the  size  of  the  beam  in 
the  crystal  produced  the  results  shown  in  Figure  4.17.  The  solid  line  reprepresents 
the  spot  size  in  the  laser  crystal  scaled  up  by  a  factor  of  1.1  over  the  value  calculated 
using  the  cavity  configuration  information.  The  curve  shows  that  increasing  the 
spot  size  decreases  the  pulse  energy  and  increases  the  pulse  length.  An  increase 
in  the  spot  size  increases  the  mode  volume  which  intersects  the  laser  crystal.  An 
increase  in  the  mode  volume  increases  the  threshold  population  inversion  (which  can 
be  seen  by  examining  Eq  (2.7)  on  page  2-6).  A  beam  size  of  1.9  times  larger  than 
the  value  calculated  from  the  cavity  geometry  produces  the  correct  pulse  lengths 
but  the  corresponding  pulse  energies  are  too  low.  Using  this  beam  size,  the  value  of 
the  ratio  was  adjusted  until  the  calculated  energy  matched  the  measured 

values.  The  value  of  the  MiisisnsiL  ratio  was  found  to  be  5.27  and  the  results  are 

“total 

plotted  in  Figure  4.18. 
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Pulse  energy(microJoules)  Pulse  width  (ns) 


Figure  4.17  Model  results  for  pulse  width  and  pulse  energy  as  functions  of  pulse 
repetition  frequency.  The  multiple  curves  demonstrate  the  effect  of 
changing  the  spot  size  in  the  crystal. 
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Pulse  energy(microJoules)  Pulse  width  (ns) 


Figure  4.18  Model  results  for  pulse  width  and  pulse  energy  as  functions  of  pulse 
repetition  frequency.  The  multiple  curves  demonstrate  the  effect  of 
changing  amount  of  absorbed  power  while  maintaining  the  same  ratio 
of  absorption  efficiency  to  total  cavity  loss.  The  value  of  use(f 

to  produce  this  plot  was  5.27  (empirical). 
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4-6  Dual  wavelength  Q-switched  results 

This  section  begins  by  outlining  the  control  sequence  for  the  dual  wavlength 
cavity.  Examples  of  the  output  which  reflect  this  control  sequence  are  presented. 
Two  cavity  lengths  were  analyzed  and  the  measurement  results  are  presented  in  a 
series  of  tables.  Most  of  the  data  for  Nd:YAG  is  spread  out  in  parameter  space 
and  the  initial  plan  was  to  use  this  data  to  optimize  the  model  which  could  then 
predict  the  optimum  delay  time,  and  cavity  configuration.  The  development  never 
reached  that  stage.  The  results  do,  however,  illustrate  some  key  concepts  which  are 
important  for  designing  this  type  of  laser  for  a  laser  radar  system. 

Figure  4.19  on  page  4-37  shows  the  timeline  for  the  control  sequence  as  well 
as  some  of  the  terms  used  to  describe  the  operation  of  the  dual  wavelength  cavity. 
When  power  is  applied  to  the  acousto-optic  Q-switches,  they  deflect  a  portion  of  the 
incident  radiation  making  the  cavity  high-loss.  Once  a  Q-switch  is  turned  “off” ,  it 
no  longer  deflects  energy  and  the  cavity  returns  to  low-loss. 

The  experimental  set  up  for  these  pulsed  measurements  is  shown  in  Figure  4.11 
on  page  4-21  and  the  procedure  for  aligning  the  cavity  is  described  in  Appendix  A. 
A  plot  of  the  shape  of  the  1.06/xm  and  1.32/im  pulses  is  shown  in  Figure  4.20  on 
page  4-38.  It  should  be  noted  that  the  time  scales  for  the  two  plots  in  the  figure  are 
different.  The  outputs  from  the  two  detectors,  showing  the  1.3fxm  and  the  1.06/im 
pulses  are  shown  in  Figure  4.21  on  page  4-39.  The  figure  shows  the  amplitude  of 
the  1.06/xm  pulse,  which  is  reflected  by  the  dichroic  beamsplitter  into  the  1.3/xm 
output  leg,  is  higher  than  the  amplitude  of  the  1.3/im  pulse.  This  highlights  one 
design  feature  that  needs  refinement  for  this  system  to  be  optimized  for  Nd:YAG. 
The  transmission  of  the  dichroic  beamsplitter  at  1.064^m  was  measured  to  be  89.5%. 
This  is  significantly  lower  than  the  design  specification  of  95%  for  unpolarized  light. 
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Figure  4.19  Timeline  showing  the  control  sequence  for  the  operation  of  the  dual 
wavelength  cavity. 
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Detector  Response  (arb.  units.)  Detector  Response  (arb.  units.) 


time  (nanoseconds) 

Figure  4.20  Pulse  shape  of  the  1.06/xm  and  1.32/im  pulses  from  the  Nd:YAG  output. 
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Detector  Response  (a.  u.) 


Nd:YAG 


Figure  4.21  Example  output  of  the  dual  wavelength  cavity.  The  top  curve  shows 
the  output  of  the  1.3/im  cavity  measured  with  the  ET3000  detector 
and  it  shows  the  presence  of  both  the  1.3/im  pulse  and  the  1.06/im 
pulse  which  is  reflected  by  the  dichroic  into  the  1.3/im  output  leg.  The 
bottom  curve  shows  the  output  of  the  1.06/im  cavity  measured  with 
the  FND-100  detector. 
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4-6.1  Energy  Crossover.  Table  4.14  on  page  4-42  contains  the  results  from 
Nd:YLF  operating  in  the  dual  wavelength  mode  of  operation.  The  pulse  energy  and 
pulse  length  data  is  plotted  as  a  function  of  the  delay  time  in  Figure  4.22  on  page  4- 
41.  As  the  delay  time  increases,  it  increases  the  population  inversion  for  the  1.047pm 
pulse  and  decreases  it  for  the  1.3pm  pulse.  This  decrease  in  the  initial  population 
inversion  decreases  the  portion  of  initial  population  which  is  extracted  according  to 
Eq  (2.22).  Population  left  over  by  the  inefficiency  of  the  1.32pm  provides  additional 
population  inversion  for  the  1.05pm  pulse  which  increases  its  extraction  efficiency. 
By  increasing  the  extraction  efficiency  of  the  1.05pm  pulse,  the  initial  population 
inversion  for  the  1.32pm  pulse  decreases  even  further,  and  so  the  whole  process 
reinforces  a  decrease  in  the  amount  of  1.32pm  energy.  In  order  to  have  a  stable 
system  which  has  equivalent  pulse  energies  at  1.32pm  and  1.05pm,  the  pumping 
rate  must  be  sufficient  at  the  given  pulse  repetition  frequency  to  provide  enough 
initial  population  inversion  relative  to  the  threshold  inversion  level  for  the  1.32pm  to 
efficiently  extract  the  population  inversion.  The  population  inversion  should  be  at 
least  twice  the  threshold  inversion  level  for  the  1.3pm  pulse.  The  relationship  of  the 
extraction  efficiency  to  the  ratio  of  initial  population  inversion  to  threshold  inversion 
level  ( Ui/nth )  in  Figure  2.5  on  page  2-12  shows  that  for  ratios  of  less  than  2.0,  the 
extraction  efficiency  starts  dropping  dramatically. 

In  designing  a  laser  radar  system,  the  energy  of  the  two  wavelengths  would  be 
adjusted  to  account  for  differences  in  atmospheric  propagation,  target  reflectivity, 
detector  responsivity,  and  other  performance  aspects  considered  in  a  system  design 
and  trade  off  analysis.  From  Figure  4.22,  it  is  clear  that  for  a  fixed  repetition  rate, 
the  amount  of  delay  between  the  1.3pm  pulse  and  the  1.064pm  pulse  can  be  used  to 
change  the  relative  amount  of  energy  in  the  two  pulses. This  can  be  seen  more  clearly 
in  Figure  4.23  where  the  pulse  energy  from  both  wavelengths  from  Figure  4.22  are 
plotted  on  the  same  graph.  The  two  plots  from  the  two  wavelenths  actually  cross. 
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1047  nm 


1320  nm 


Figure  4.22  Dual  wavelength  operation  of  Nd:YLF.  The  repetition  rate  for  each 
wavelength  was  fixed  at  5.0kHz  and  the  absorbed  pump  power  was  950 
mW. 
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Table  4.14  Q-switched  Nd:YLF  dual  wavelength  measurements-  short  cavity. 


Pump 

Power 

Inc. 

(W) 

Pump 

Power 

Abs. 

(W) 

Pulse 

Rep. 

Rate 

(Hz) 

Delay 

time 

(/xs) 

1.06/xm 

Pulse 

Energy 

0<J) 

1.06/xm 

Pulse 

Width 

(ns) 

1.34/xm 

Pulse 

Energy 

0*J) 

1.34/xm 

Pulse 

Width 

(ns) 

0.965 

0.92±0.03 

5000 

85 

11.2±0.9 

109±7 

28±1 

220±10 

0.965 

0.92±0.03 

5000 

95 

12.8±0.9 

109±3 

27±1 

237±7 

0.965 

0.92±0.03 

5000 

105 

14.8±1 

89±8 

26±1 

240±10 

0.965 

0.92±0.03 

5000 

115 

19.5±0.8 

69±4 

22.9±0.9 

280T10 

0.965 

0.92±0.03 

5000 

125 

23.1±1 

56±7 

20±1 

269±7 

0.965 

0.92±0.03 

5000 

135 

25.5±1 

55±3 

18.6T0.9 

320±20 

Changing  the  delay  time  can  compensate  for  the  differences  in  stimulated  emis¬ 
sion  cross  section  as  well  as  losses  between  the  two  cavities.  The  ability  to  adjust  the 
pulse  energies  using  the  delay  time  validates  the  use  of  the  stimulated  emission  cross 
section  -  fluorescence  lifetime  product  as  a  metric  for  choosing  the  optimum  material 
for  this  dual  wavelength  laser.  Any  defficiency  in  the  performance  of  one  transition 
can  be  made  up  by  performance  in  the  other  simply  by  changing  the  relative  amount 
of  time  the  population  inversion  is  allowed  to  build  up  for  each  pulse.  Different 
materials  with  comparable  values  of  this  metric  would  have  to  be  evaluated  for  their 
ability  to  achieve  higher  repetition  rates.  The  energy  in  the  materials  with  the  longer 
lifetimes  will  begin  to  drop  sooner  as  the  repetition  rate  is  increased  (3:623) 

4-6.2  Cavity  length  comparison.  Measurements  were  made  at  with  the 
Nd:YAG  crystal  at  two  different  cavity  lengths.  The  longer  cavity  lengths  were  on 
the  order  of  18cm  long  and  the  shorter  cavity  lengths  were  on  the  order  of  13cm 
long.  The  length  of  the  short  cavity  was  determined  by  reducing  the  distances 
between  components  as  much  as  possible  within  the  constraints  of  the  hardware.  The 
separation  of  the  different  components  in  the  short  cavity  are  shown  in  Table  4.16 
(see  Figure  4.11  on  page  4-21  for  the  definitions)  and  the  laser  measurements  are 
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Pulse  Energy  (microJoules) 


Figure  4.23  Dual  wavelength  operation  of  Nd:YLF  showing  energy  crossover.  The 
repetition  rate  for  each  wavelength  was  fixed  at  5.0kHz  and  the  ab¬ 
sorbed  pump  power  was  920  mW. 


Table  4.15  Q-switched  Nd:YAG  dual  wavelength  measurements-short  cavity. 


Pump 

Power 

Inc. 

(W) 

Pump 

Power 

Abs. 

(W) 

Pulse 

Rep. 

Rate 

(Hz) 

Delay 

time 

(fis) 

1.06/xm 

Pulse 

Energy 

0J) 

1.06/nn 

Pulse 

Width 

(ns) 

1.34/xm 

Pulse 

Energy 

w 

1.34/xm 

Pulse 

Width 

(ns) 

1 

0.83±0.03 

2000 

100 

33±1 

31.5±0.4 

26.7±0.8 

183±1 

1 

0.83±0.03 

2000 

150 

40±1 

32.9±0.7 

23.2±0.7 

199±4 

0.83±0.03 

2000 

225 

43±1 

28.2±0.6 

19.3±0.6 

274±5 

0.83±0.03 

2880 

100 

32±1 

32.2±0.3 

16.3±0.5 

338±3 

Table  4.16  Cavity  configuration  for  the  dual- wavelength  Q-switched  Nd:YAG  data 
-short  cavity. 


li=1.4cm 

14=0. 7cm 

l2=3.5cm 

15=2. 4cm 

13=3. 5cm 

l6=1.9cm 

Output  coupler 

Cavity  length 

1.32/xm 

T=4%,  20cmcc 

11.4cm 

1.06/rm 

T=5%,  20cmcc 

12.2cm 

shown  in  Table  4.15.  The  separations  of  the  components  in  the  long  cavity  are  shown 
in  Table  4.10  back  on  page  4-23  and  the  laser  measurements  are  shown  in  Table  4.17. 


The  effect  of  the  changing  the  cavity  length  is  shown  in  Figure  4.24.  Shortening 
the  cavity  length  from  18cm  to  13cm  has  a  dramatic  impact  on  the  pulse  length.  The 
pulse  lengths  for  both  wavelengths  are  reduced  by  a  factor  of  two.  This  supports  the 
conclusion  that  the  shortest  possible  cavity  is  the  optimum  one  for  producing  short 
pulses(3:622).  Shortening  the  cavity  length  decreases  the  round  trip  time,  and  as  a 
result,  the  photon  lifetime  is  decreased  also.  When  this  happens,  the  photons  in  the 
cavity  provide  their  feedback  to  the  gain  medium  faster  causing  the  pulse  to  build  and 


4-44 


Table  4.17  Q-switched  Nd:YAG  dual  wavelength  measurements-long  cavity. 


Pump 

Power 

Inc. 

(W) 

Pump 

Power 

Abs. 

(W) 

Pulse 

Rep. 

Rate 

(Hz) 

Delay 

time 

(ps) 

1.06/on 

Pulse 

Energy 

M 

1.06/zm 

Pulse 

Width 

(ns) 

1.34/xm 

Pulse 

Energy 

(/*J) 

1.34/zm 

Pulse 

Width 

(ns) 

1.45 

0.92  ±  0.03 

500 

100 

42±4 

49.8±0.9 

40±2 

211±2 

1.45 

0.92  ±  0.03 

500 

150 

44±2 

43.4  ±0.3 

45±1 

216±2 

1.295 

0.823±  0.03 

500 

200 

47±2 

45.9  ±0.5 

35±1 

288±7 

1.3 

0.83±  0.03 

1000 

200 

45±2 

46.8±0.5 

33  ±1 

301±2 

lEB 

0.83  ±  0.03 

1000 

150 

42±1 

64.0±0.4 

34±1 

267±7 

mm 

0.83±  0.03 

2000 

100 

38  ±2 

52.1±0.5 

21.7±0.7 

410±10 

0.83±  0.03 

2000 

150 

38  ±2 

48.1±0.6 

19.2±0.6 

504±4 

EB 

0.76  ±  0.03 

2000 

100 

36  ±1 

54  ±1 

10.5  ±0.3 

790±20 

sweep  out  the  gain  region  faster  than  in  a  longer  cavity.  This  increased  photon  flux 
increases  the  efficiency  of  extracting  the  stored  energy  through  stimulated  emission. 

4-6.3  Material  comparison.  A  set  of  measurements  was  also  made  with 
Nd:YLF  in  a  short  cavity.  The  cavity  lengths  for  Nd:YAG  (Table  4.16)  were  ap¬ 
proximately  1cm  shorter  than  they  were  in  Nd:YLF  (Table  4.18)  because  the  angle 
of  the  dichroic  beamsplitter  was  different  for  the  two  materials  causing  the  relative 
position  of  the  Q-switch  mounts  to  change.  The  laser  measurements  for  Nd:YLF 
are  shown  in  Table  4.19.  The  measurement  results  for  Nd:YAG  were  shown  in  the 
previous  subsection  in  Table  4.15. 

The  Nd:YAG  and  Nd:YLF  measurements  made  with  the  short  cavity  were  to 
provide  a  direct  comparison  of  the  two  materials  using  the  shortest  possible  cavity. 
Unfortunately,  there  are  several  differences  between  the  results  from  the  two  ma¬ 
terials  which  make  a  direct  comparison  impossible.  The  analysis  in  section  4.4  for 
the  1.06/zm  cavity  makes  the  calculation  for  the  absorbed  pump  power  uncertain, 
it  indicates  that  the  level  of  absorbed  power  for  Nd:YAG  might  be  greater  than 
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Pulse  Energy  (microJoules)  Pulse  length  (ns) 


Figure  4.24  Dual  wavelength  operation  of  Nd:YAG  comparing  the  results  from  two 
different  cavity  lengths.  The  plots  on  the  left  hand  side  are  for  the 
1.06//m  cavity  and  the  plots  on  the  right  hand  side  are  for  the  1.32/mi 
cavity.  The  pulse  repetition  frequency  (for  both  wavelenths)  is  fixed  at 
2000Hz. 


Table  4.18  Cavity  configuration  for  the  dual-wavelength  Q-switched  Nd:YLF  data 
-short  cavity. 


li=1.4cm 

14=0. 7cm 

12=4. 1cm 

15=2. 1cm 

13=4. 6cm 

16=2. 5cm 

Output  coupler 

Cavity  length 

1.32/im 

T=4%,  20cmcc 

13.1cm 

1.06//m 

T=30%,  20cmcc 

13.0cm 
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Table  4.19  Q-switched  Nd:YLF  dual  wavelength  measurements-short  cavity. 


Power 

Inc. 

(W) 

Power 

Abs. 

(W) 

Rep. 

Rate 

(Hz) 

Delay 

time 

(/xs) 

1.047/xm 

Pulse 

Energy 

(ix3) 

1.047/xm 

Pulse 

Width 

(ns) 

1.32/xm 

Pulse 

Energy 

0*j) 

1.32/xm 

Pulse 

Width 

(ns) 

0.959 

0.76±0.02 

2000 

100 

34±1 

51±2 

12.3  ±0.5 

460±10 

0.959 

0.76±0.02 

2000 

150 

39±1 

44.4±0.5 

9.9±0.4 

666±7 

0.959 

0.76±0.02 

2000 

225 

41±2 

44.8±0.9 

9.8±0.4 

1020±40 

1.046 

0.82±0.03 

2000 

100 

40±2 

44±1 

15.3±0.6 

410±10 

1.046 

0.82±0.03 

2000 

150 

48±2 

45  ±1 

13.8±0.5 

470±20 

1.046 

0.82±0.03 

2000 

225 

52±2 

38.4±0.4 

12.7±0.7 

560±20 

originally  calculated.  This  means  that  the  two  materials  were  pumped  at  different 
rates.  Not  only  was  the  pumping  rate  different,  but  the  output  coupler  values  were 
different.  The  gain  in  the  Nd:YLF  1.047/xm  laser  was  so  high  that  the  Q-switch 
could  not  hold  off  lasing  between  pulses.  The  output  coupler  reflectivity  was  then 
reduced  from  95%  to  70%  where  the  1.06/xm  Q-switch  could  hold  off  the  cavity  and 
this  is  the  point  where  the  data  was  taken.  When  the  output  coupler  reflectivity  for 
the  1.06/xm  cavity  was  reduced  from  95%  to  70%  in  Nd:YAG,  the  loss  was  too  high 
and  the  cavity  would  not  lase  at  that  pulse  repetition  rate  and  pump  power.  As  a 
result,  the  reflectivity  of  output  coupler  for  the  Nd:YLF  data  was  70%,  whereas  the 
reflectivity  of  the  output  coupler  for  the  Nd:YAG  data  was  95%.  The  results  from 
the  two  sets  of  measurements  are  plotted  in  Figure  4.25. 

Figure  4.25  again  demonstrates  the  interdependence  of  the  two  wavelengths.  As 
the  amount  of  time  the  population  is  allowed  to  build  for  the  1.06/xm  pulse  (the  delay 
time)  is  increased,  the  amount  of  population  inversion  increases  for  the  1.06/xm  pulse 
at  the  expense  of  the  1.3/xm  pulse.  The  reason  for  the  pulse  length  of  the  Nd:YLF 
4F3/2  -  4Iii/2  transition  being  longer  than  the  pulse  length  for  Nd:YAG  when  Nd:YLF 
has  a  greater  pulse  energy  is  not  immediately  obvious.  One  possibility  is  that  the 
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Pulse  Energy  (microJoules)  Pulse  length  (ns) 


Figure  4.25  Dual  wavelength  operation  of  Nd:YLF  compared  with  Nd:YAG.  The 
plots  on  the  left  hand  side  are  for  the  1.06//m  cavity  and  the  plots  on 
the  right  hand  side  are  for  the  1.32/im  cavity.  The  pulse  repetition 
frequency  (for  both  wavelenths)  is  fixed  at  2.0kHz. 
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1.06/jm  Q-switch  did  not  have  enough  hold  off  to  keep  that  cavity  from  lasing  between 
pulses.  During  the  measurement  the  Q-switch  power  was  at  its  absolute  maximum, 
and  it  was  not  clear  from  the  amplitude  of  the  1.047/im  Q- switched  pulse  that  the 
Q-switch  was  able  to  fully  hold  off  the  gain  in  the  cavity.  The  Q-switch  was  designed 
for  unpolarized  light.  The  relative  efficiency  of  the  device  for  horizontally  polarized 
energy  was  never  evaluated. 

4.7  Error  analysis 

4-7.1  Uncertainty  in  the  pump  power  absorbed  in  the  laser  crystal.  The 
uncertainty  in  the  amount  of  absorbed  power  results  from  the  uncertainty  in  the 
absolute  accuracy  of  the  power  meter;  the  fluctuation  in  the  meter  reading  (which 
was  recorded  for  each  measurement);  the  uncertainty  in  the  transmission  of  the  high 
reflector;  and  the  uncertainty  in  the  fraction  of  incident  energy  actually  absorbed 
by  the  laser  crystal.  The  setup  for  this  measurement  is  shown  in  Figure  4.3.  The 
expression  used  to  calculate  the  aborbed  pump  power  from  the  incident  pump  power 
Pine  was 


PAbs  =  VHRVAbsVBPinc  (4-8) 

The  values  for  ?}hr  and  rjAbs  were  taken  from  Tables  4.1  and  4.3  respectively. 
The  value  for  t]b  is  assumed  to  be  1,  the  pump  beam  was  focused  to  a  spot  size 
smaller  than  the  spot  size  of  the  laser  mode  inside  the  laser  crystal. 

The  accuracy  of  the  power  meters  <Jaccu  used  to  make  the  power  measurements 
was  ±3%.  The  meter  reading  for  Pmddent  typically  fluctuated  ±5mW  during  a  mea¬ 
surement.  This  fluctuation  is  included  in  the  uncertainty  calculation.  The  equation 
used  to  calculate  the  uncertainty  in  the  power  absorbed  into  the  laser  crystal  <JpAba 


4-49 


was 


where  aacc  is  the  absolute  accuracy  of  the  power  meter,  (Tfiuct  is  the  fluctuation  in 
the  meter  reading  ,  <JtHr  is  the  uncertainty  in  the  transmission  of  the  high  reflector, 
and  <JVAbs  is  the  uncertainty  in  the  fraction  of  power  absorbed  by  the  crystal. 

Representative  numbers  for  the  fractions  inside  the  expression  in  Eq  (4.9)  for 
Nd:YAG  pumped  with  1.0  Watt  from  the  Ti:Sapphire  at  808.5  nm  are  given  below 


a 


decline 

~p. 

1  me 

(tfiuct 


p. 

1  me 
aTHR 


Thr 

°VAbB 


VAbs 


0.03 

0.005 

0.0135 

0.011 


A  comparison  of  these  values  indicates  that  the  first  one,  from  the  absolute  accuracy 
of  the  power  meter,  is  by  a  factor  of  3  the  greatest  source  of  uncertainty  in  the 
calculation  of  the  amount  of  pump  power  absorbed  in  the  laser  crystal. 


7.2  Uncertainty  calculation  for  the  laser  pulse  energy.  The  equation  used 
to  calculate  the  pulse  energy  from  the  average  power  measured  with  the  Scientech 
model  AD30  power  meter  was 


Wout  = 


avg 


TlpU  foldF  fold  (1-T  0C)PRF 


(4.10) 


where  Pavg  is  the  average  power  measured  with  the  power  meter,  TLp  is  the  trans¬ 
mission  of  the  long  pass  filter  at  the  laser  wavelength,  T jm  is  the  reflectivity  of 
the  protected  silver  fold  mirrors,  Toe  is  the  reflection  off  of  the  back  surface  of  the 
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output  coupler  (the  output  couplers  were  not  coated  on  the  back  side),  and  PRF  is 
the  pulse  repetition  frequency. 


The  uncertainty  in  the  pulse  energy  measurements  results  from  the  uncertainty 
in  the  absolute  accuracy  of  the  power  meter;  the  fluctuation  in  the  meter  reading 
(which  was  recorded  for  each  measurement);  the  uncertainty  in  the  transmission  of 
the  long  pass  filter  used  to  block  the  residual  pump  beam;  and  the  uncertainty  in 
the  reflectivity  of  the  protected  silver  fold  mirror(s)  used  to  direct  the  beam  into  the 
power  meter.  The  reflectivity  off  of  the  back  surface  of  the  optic  was  not  measured 
but  assumed  to  be  a  constant  5%  with  no  associated  uncertainty.  The  uncertainty 
in  the  pulse  repetition  frequency  was  also  disregarded  since  it  was  less  than  one  part 
in  one  thousand,  an  order  of  magnitude  less  than  the  other  sources  of  error.  The 
equation  used  to  calculate  the  uncertainty  in  the  pulse  energy  was 


where  Pavg  is  the  average  power  of  the  pulsed  laser  measured  with  the  power  meter. 

Representative  calculations  for  Nd:YAG  at  1.06/im  associated  with  the  top 
entry  in  Table  4.8  on  page  4-22  are  shown  below 


® accP ivg 


avg 


G fluct 

~P 

1  avg 

fold 
r  fold, 
aTLP 

Tlp 


0.03 

0.02 

0.0052 

0.0094 


Comparing  these  values  indicates  that  the  first  value,  from  the  absolute  accuracy  of 
the  power  meter,  is  the  greatest  source  of  uncertainty  in  the  calculation  of  the  amount 
of  pump  power  absorbed  in  the  laser  crystal.  When  the  average  power  produced  by 
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the  laser  was  less  than  20mW,  the  fluctuation  in  the  meter  reading  started  having  a 
significant  impact  on  the  uncertainty  values.  This  generally  occured  at  low  repetition 
rates  or  low  pump  powers. 


4-52 


V.  Conclusions  &  Recommendations 

5. 1  Conclusions 

A  figure  of  merit  for  choosing  the  best  laser  crystal  from  the  large  number  of 
new  Nd-based  materials,  is  the  sum  of  the  stimulated  emission-fluorescence  lifetime 
product  for  the  two  transitions  that  produce  the  two  wavelengths.  Due  to  both  the 
large  variability  and  the  gaps  in  the  published  values  for  these  material  parame¬ 
ters,  the  best  method  for  choosing  the  optimum  material  is  to  evaluate  their  relative 
performance  in  an  identical  laser  setup.  Based  on  the  published  information  avail¬ 
able,  Nd:SFAP  would  be  the  optimum  material  for  this  application  because  the  sum 
of  the  stimulated  emission  cross  section-fluorescence  lifetime  product  from  the  two 
wavelengths  is  23.2  x  10~17/is  cm2,  which  is  higher  than  any  other  material  surveyed. 

The  absorption  spectrum,  fluorescence  lifetime,  and  emission  spectrum  of  Nd:YAG, 
Nd:SVAP,  and  Nd:YLF  were  measured.  The  peak  absorption  in  the  band  of  useful 
diode  wavelengths  between  780nm  and  815nm  for  both  Nd:YAG  and  Nd:SVAP  was 
approximately  808. 5nm,  whereas  the  strongest  absorption  for  Nd:YLF  occurred  at 
792. lnm.  The  measured  fluorescence  lifetimes  for  these  materials  was  237/xs,  216/zs 
and  504/7, s  respectively,  and  all  are  within  the  range  of  published  values. 

A  dual  wavelength  Q-switched  laser  cavity  has  been  successfully  designed  and 
assembled  to  evaluate  the  different  possible  laser  materials.  The  cavity  design  was 
based  on  a  “Y”  cavity  designed  by  Tomaschke  (33:  252)  to  achieve  lasing  at  two 
wavelengths  simultaneously.  Initial  characterization  has  been  achieved  for  Nd:YAG 
and  Nd:YLF.  The  results  indicate  that  for  a  fixed  pulse  repetition  frequency,  the 
delay  time  can  be  used  to  adjust  the  relative  energy  between  the  1.06pm  and  1.32pm 
pulses.  This  ability  validates  using  the  sum  of  the  fluorescence  lifetime-stimulated 
emission  cross  section  product  from  the  two  transitions  as  a  metric  for  choosing  the 
best  material.  Any  deficiency  in  the  performance  of  one  transition  can  be  made  up 
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by  performance  in  the  other  simply  by  changing  the  relative  amount  of  time  the 
population  inversion  is  allowed  to  build  up  for  each  pulse. 

The  shortest  pulses  were  achieved  with  both  Nd:YAG  and  Nd:YLF  using  the 
shortest  physical  cavity  length.  When  the  cavity  length  was  decreased  from  18cm  to 
13cm,  the  pulse  length  for  each  wavelength  decreased  by  a  factor  of  2. 

Of  the  two  materials  evaluated,  Nd:YLF  was  clearly  the  superior  material.  The 
gain  was  so  high  with  this  material  that  the  reflectivity  of  the  output  coupler  had  to 
be  decreased  from  95%  to  70%  just  so  that  the  1.06/im  Q-switch  could  hold  off  the 
cavity. 

The  model  developed  using  the  standard  differential  equations  was  fit  to  the 
Nd:YAG  data  at  1.06/un.  The  modifications  which  were  required  to  make  the  model 
fit  the  data  indicate  the  presence  of  thermal  lensing  in  this  material.  An  improved 
model  would  include  the  effects  of  thermal  lensing  as  a  function  of  pulse  repetition 
rate  and  energy  storage  in  the  material.  The  importance  of  thermal  lensing  was  also 
seen  in  the  single  wavelength  results  for  Nd:YLF.  The  pulse  energy  actually  started 
to  decrease  with  decreasing  repetition  rate  instead  of  continuing  to  increase  as  would 
normally  be  exptected. 

5.2  Recommendations 

The  first  and  most  obvious  recommendation  is  to  obtain  a  sample  of  NdrSFAP 
for  evaluation.  Based  on  the  published  spectroscopic  information,  this  material 
should  be  the  best  out  of  those  surveyed  for  this  two  wavelength  application.  Based 
on  the  information  in  Table  3.8  on  page  3-5,  Nd:YLF  is  the  next  best  material  to 
Nd:SFAP.  The  performance  of  these  materials  in  terms  of  pulse  energy  and  pulse 
width  should  be  evaluated  up  to  the  pulse  repetition  rates  necessary  to  achieve  the 
needed  frame  rate  for  a  dual  wavelength  imaging  laser  radar.  In  addition,  these 
materials  should  be  evaluated  not  only  for  their  pulse  energy  and  pulse  width  but 
for  their  absorption  width  at  the  laser  diode  pump  wavelength  (to  evaluate  their 
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suitability  for  laser  diode  pumping),  thermal  stability,  and  beam  quality  as  well. 
All  of  these  parameters  have  an  impact  on  the  performance  of  a  laser  radar  system 
constructed  from  these  materials. 

The  cavity  element  which  caused  the  greatest  loss  and  proved  to  be  the  greatest 
detriment  to  operation  on  the  4F3/2  -  4In/2  transition  was  the  dichroic  beamsplitter. 
This  element  had  a  transmission  of  90%  in  the  Nd:YAG  cavity.  One  possible  solution 
is  to  purchase  a  higher  quality  optic  with  a  better  design.  Another  solution  is  to 
modify  the  cavity  design  and  use  a  polarization  beamsplitter. 

Since  both  of  the  laser  transitions  in  all  of  the  laser  crystals  have  a  significant 
stimulated  emission  cross  section  in  both  polarization  directions,  polarizing  elements 
could  be  placed  in  each  of  the  different  legs  of  the  cavity.  The  1.3/xm  polarizer  could 
be  oriented  to  allow  only  the  vertical  to  pass.  This  vertical  polarization  would  be 
S-polarized  upon  incidence  the  dichroic  beamsplitter  and  so  its  reflection  would  be 
maximized.  The  polarizer  in  the  1.06/im  cavity  could  be  oriented  to  allow  only 
horizontal  polarization  to  pass.  This  horizontal  polarization  would  be  P-polarized 
upon  incidence  to  the  dichroic  beamsplitter  which  would  maximize  its  transmis¬ 
sion.  If  both  legs  then  contain  optics  to  force  opposite  polarizations,  the  dichroic 
beamsplitter  could  be  replaced  with  a  polarization  beamsplitter.  This  polarization 
beamsplitter  would  have  to  be  designed  to  operate  at  both  wavelenths  in  the  near 
IR,  but  a  beamsplitter  made  from  calcite  would  give  an  extinction  ratio  between  the 
two  polarizations  of  100,000  to  1.  The  potential  increase  in  performance  would  have 
to  be  evaluated  in  terms  of  the  additional  optical  loss  caused  by  having  an  additional 
intracavity  element  as  well  as  the  added  complexity. 

One  characteristic  of  the  output  which  was  not  examined  here  was  the  wave¬ 
length  content  of  the  Q-switched  output.  The  CW  measurements  made  with  Nd:YAG 
and  Nd:YLF  indicate  the  presence  of  more  than  one  wavelength  at  1.3/rai.  The  tar¬ 
get  discrimination  capability  of  a  dual  wavelength  laser  radar  system  depends  not 
only  on  the  wavelength  of  each  output  but  also  on  the  separation. 
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Temperature  effects  on  the  performance  of  the  system  need  to  be  studied. 
Other  researchers  have  found  that  reducing  the  temperature  of  certain  crystals  helps 
increase  the  gain  of  the  laser  system  allowing  higher  energies  and  shorter  pulses 
(22:886).  The  design  of  the  TE  cooled  laser  crystal  mount  incorporates  enough 
capacity  to  make  measurements  down  to  at  least  0°  C.  In  order  to  make  these 
measurements,  an  enclosure  over  the  laser  crystal  mount  would  have  to  be  fabricated 
and  connected  to  a  dry  nitrogen  purge  in  order  to  keep  moisture  from  condensing 
on  the  mount  and  crystal. 

The  system  would  benefit  greatly  from  an  enclosure  which  could  be  kept  in 
place  during  operation.  The  spot  size  of  the  laser  beam  in  the  laser  crystal  is  on 
the  order  of  150/mi.  This  is  on  the  same  order  magnitude  of  some  of  the  larger  dust 
particles  in  the  laboratory.  Significant  variations  in  the  performance  of  the  system 
were  experienced  when  the  optics  were  not  clean  enough  and  significant  blocks  of 
time  were  spent  optimizing  the  performance  due  to  this.  To  minimize  this  impact, 
an  enclosure  which  is  purged  with  dry  nitrogen  is  necessary. 
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Appendix  A.  Laser  system  set  up 


A.l  Setup  for  the  Ti:Sapphire  Pump  source 

A.  1.1  Ti:Sapphire  laser.  The  experimental  setup  is  shown  in  Figure  4.1, 
on  page  4-2.  A  pellicle  beam  splitter  was  used  to  pull  off  a  small  portion  of  the 
pump  beam  and  send  it  to  a  Burleigh  Wavemeter  to  measure  the  optical  wavelength 
of  the  Ti:Sapphire  laser.  This  setup  proved  extremely  convenient  for  tuning  the 
Ti:Sapphire  to  the  maximum  absorbtion  of  each  material. 

The  only  purpose  of  the  two  fold  mirrors  was  to  bring  the  pump  beam  down 
to  the  elevation  of  the  laser  crystal.  The  alignment  of  the  pump  beam  to  the  laser 
crystal  was  performed  without  the  lens  in  place  using  the  second  fold  mirror.  The 
30cm  focal  length  lens  used  to  focus  the  Ti:Sapphire  pump  beam  down  into  the 
crystal  had  a  tendency  to  steer  and  distort  the  pump  beam  if  it  was  not  aligned 
normal  to  the  beam.  After  using  the  second  fold  mirror  to  align  the  beam  to  the 
center  of  the  crystal,  the  lens  was  placed  in  the  beam  path  and  positioned  so  that 
the  laser  spot  struck  the  center  of  the  laser  crystal.  It  was  then  aligned  normal  to  the 
pump  beam  by  making  sure  the  back  reflections  off  of  the  lens  surfaces  were  headed 
back  into  the  aperture  of  the  Ti:Sapphire. 

The  transmission  of  the  high  reflector  cavity  mirror,  at  the  pump  wavelength, 
proved  to  be  lower  than  expected.  Although  the  back  surface  (the  surface  the  pump 
beam  strikes  first)  is  anti-reflective  coated  at  the  pump  wavelength,  the  second  sur¬ 
face,  the  dielectric  high  reflector,  reflects  approximately  25%  of  the  incident  pump. 
When  this  surface  was  normal  to  the  pump  beam — which  it  was  in  most  cases  for  the 
most  efficient  use  of  the  pump  beam  in  the  laser  mode  volume — this  power  would 
reflect  back  into  the  Ti:Sapphire  cavity  and  cause  the  output  of  the  Ti:Sapphire 
to  go  into  a  TEMio  mode  and  would  also  make  the  output  power  fluctuate.  The 
neodymium  laser  output  power  would  then  start  to  fluctuate,  sometimes  by  more 
than  50%.  A  more  robust  optical  design  would  include  both  a  high  reflector  which  is 
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more  transmissive  in  the  region  around  808nm  and  an  optical  isolator  in  the  pump 
beam  path. 

The  neutral  density  filter  wheel  has  a  continuous  neutral  density  gradient  and 
it  was  used  to  vary  the  Ti:Sapphire  pump  power.  The  meter  used  to  measure  the 
power  was  a  Scientech  model  365  with  a  one  inch  disc  calorimeter. 

The  pump  source  for  the  dual-wavelenth  laser  is  an  Argon-pumped  Ti:Sapphire 
laser.  The  wavelength  was  fully  tunable  over  a  much  wider  range  than  the  792nm 
to  808. 5nm  range  used  in  these  experiments.  The  power  did  not  vary  appreciably 
while  the  wavelength  was  being  tuned  over  this  range  making  it  ideal  for  tuning  to 
the  absorbtion  peak  of  each  material. 

The  waist  size  and  location  of  the  Ti:Sapphire  were  measured  using  the  Co¬ 
herent  Modemaster.  The  results  showed  a  beam  waist  located  10±  0.5cm  behind 
the  output  port  of  the  laser  with  a  beam  waist  diameter  of  720±10/xm.  This  mea¬ 
surement  was  used  to  calculate  the  placement  of  the  30cm  focal  length  lens  used  to 
focus  the  pump  beam  down  into  the  crystal.  The  distance  between  the  Ti:Sapphire 
source  beam  waist  and  the  laser  crystal  was  fixed  at  1.393±0.002m.  A  translation 
matrix  was  set  up  in  the  mathematical  package  Mathcad  which  had  variables  for 
translation  from  the  Ti:  Sapphire  through  the  lens  and  the  pump  input  mirror  which 
was  a  40  centimeters  concave  dielectric-coated  mirror.  The  translation  matrix  had 
variables  for  the  distance  from  the  waist  of  the  Ti: Sapphire  laser  to  the  lens,  and 
for  the  distance  from  the  lens  to  the  laser  crystal.  The  software  then  iterated  until 
it  found  a  location  of  the  lens  which  produced  a  147/xm  spot  size  in  the  laser  crys¬ 
tal.  This  location  was  41.5cm  from  the  input  face  of  the  laser  crystal.  The  value  of 
147/xm  was  chosen  in  order  to  match  the  calculated  spot  size  of  the  laser  mode  in 
the  center  of  the  laser  crystal.  After  the  lens  was  physically  placed  in  the  calculated 
location,  the  beam  size  was  again  checked  with  the  Coherent  Modemaster  and  found 
to  be  146±5/xm. 
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A.  1.2  Laser  crystal  mount.  The  mounts  that  held  the  laser  crystal  samples 
were  machined  from  solid  brass.  These  mounts  were  separated  from  the  solid  copper 
base  by  three  thermo-electric(TE)  coolers.  The  coolers  had  a  sum  total  of  approx¬ 
imately  12  Watts  of  cooling  capacity  at  room  temperature.  During  the  design,  the 
most  heat  the  TE  coolers  were  envisioned  having  to  dissipate,  was  the  full  1  Watt  of 
pump  power  from  the  Ti:Sapphire  laser.  The  reason  for  the  designed-in  over  capacity 
was  two-fold.  First,  it  provides  very  consistent,  worry-free  operation;  and  second,  it 
provides  the  capability  to  do  low  (0°C)  temperature  measurements. 

A. 2  Alignment  of  the  dual  wavelength  laser  cavity 

The  first  step  in  setting  up  the  dual  cavity  was  to  clean  all  the  optics.  This  had 
to  be  repeated  later  in  the  alignment  process  depending  on  the  laser  performance. 
Dust  particles  on  the  face  of  the  laser  crystal  showed  up  distinctly  in  the  laser  output, 
causing  the  beam  to  spread  and  artifacts  to  appear  outside  the  main  beam  pattern. 

The  first  step  in  the  alignment  procedure  for  the  dual  cavity  was  to  align  the 
1.06/v.m  cavity  without  the  dichroic  or  the  Q-switch  in  the  cavity.  The  laser  output 
was  observed  on  a  screen  at  the  same  time  as  the  spot  due  to  the  residual  pump 
beam.  The  cavity  was  “walked”  so  that  the  laser  output  overlapped  the  residual 
pump  beam.  When  this  occured,  the  pump  beam’s  back-reflection  off  of  the  highly 
reflective  surface  of  the  high  reflector  was  directed  straight  back  into  the  cavity  of 
the  Ti:Sapphire,  and  the  transverse  mode  of  the  Ti:Sapphire  switched  to  TEMio-  To 
correct  this,  the  angle  of  the  high  reflector  was  then  adjusted  to  the  point  where  the 
Ti:Sapphire  switched  back  to  TEMqo  mode,  and  the  output  coupler  was  adjusted  to 
achieve  maximum  1.06/zm  power. 

The  next  step  was  to  place  the  dichroic  in  the  cavity  and  adjust  the  angle  to 
achieve  maximum  output  at  1.06/im.  The  angle  was  different  between  Nd:YAG  and 
Nd:YLF,  this  difference  most  likely  being  caused  by  the  difference  in  wavelength 
between  the  two  lasers.  The  process  of  finding  the  optimum  angle  was  an  iterative 
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one  consisting  of  finding  the  best  angle  and  then  adjusting  the  angle  of  the  output 
coupler. 

The  next  step  in  the  process  was  to  align  the  1.3/xm  cavity  with  the  dichroic  but 
not  the  Q-switch  in  the  cavity.  This  could  be  done  fairly  easily  by  letting  the  1.06/mi 
cavity  lase  while  adjusting  the  1.3/rni  output  coupler.  Alignment  was  achieved  when 
the  reflection  of  the  1.06/xm  beam  off  of  the  1.3/im  output  coupler  overlapped  the 
1.06/im  beam  reflected  from  the  back  surface  of  the  dichroic  beamsplitter. 

The  final  step  was  to  place  the  Q-switch  in  each  cavity  aligned  so  that  the 
unabsorbed  pump  beam  that  passed  through  the  laser  crystal  was  centered  on  the 
input  aperture.  The  Q-switches  were  mounted  on  4-axis  alignment  stages  so  while 
the  position  of  the  front  aperture  remained  fixed,  the  back  side,  away  from  the  laser 
crystal,  could  be  adjusted  to  achieve  maximum  power. 
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Appendix  B.  Model 


This  model  was  developed  to  configure  the  design  of  the  dual  wavelength  laser. 
The  primary  goal  in  choosing  the  radius  of  curvatures  for  the  output  couplers  was 
to  match  beam  diameter  of  the  pump  beam  and  both  laser  modes  inside  the  laser 
cavity.  Using  this  constraint,  the  output  coupler  radii  of  curvature  were  calculated 
from  a  nominal  cavity  length  of  approximately  20cm.  Once  the  nominal  size  of  the 
mode  volume  had  been  determined  by  this  geometry,  the  reflectivity  of  the  output 
couplers  were  chosen  to  achieve  at  least  twice  the  threshold  population  inversion  for 
a  constant  absorbed  pump  power  of  1  Watt. 

The  first  step  in  the  model  is  to  define  the  terms.  The  following  table  provides 
a  comprehensive  listing  of  the  variables  and  constants. 
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ll  = 

distance  from  the  high  reflector  to  the  face  of  the  laser  crystal 

h  = 

distance  from  the  laser  crystal  to  the  dichroic  beamsplitter 

I3  = 

distance  from  the  dichroic  beamsplitter  to  the  1.3/wn  Q-switch 

14  = 

distance  from  the  1.3/im  Q-switch  to  the  1.3/rni  output  coupler 

15  = 

distance  from  the  dichroic  beamsplitter  to  the  1.06/mi  Q-switch 

16  = 

distance  from  the  1.06/xm  Q-switch  to  thel.06/xm  output  coupler 

15  = 

length  of  the  gain  medium 

13  = 

length  of  the  Q-switch  (1.06/un  &  1.3/im  Q-switch  were  the  same  length) 

ns  = 

index  of  refraction  of  the  gain  medium 

n<i  = 

index  of  refraction  of  the  dichroic 

n?  = 

index  of  refraction  of  the  Q-switch 

ld  = 

distance  the  1.06/im  beam  travels  inside  the  dichroic  beamsplitter 

td  = 

thickness  of  the  dichroic  beamsplitter 

= 

The  radius  of  curvature  of  the  output  coupler 

■Rff#  = 

The  radius  of  curvature  of  the  cavity  high  reflector 

Q9 

Power  reflectivity  of  the  output  coupler 

Ti  = 

Power  reflectivity  of  the  cavity  high  reflector 
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h  = 

6.626  x  10-34  Js  Planck’s  constant 

Al.06 

Laser  wavelength  in  the  vicinity  of  1.06/rm 

Al.32 

Laser  wavelength  in  the  vicinity  of  1.32/rm 

C  — 

3xl010— 

S 

&SE  = 

Stimulated  emission  cross  section  at  the  laser  wavelength 

Tf  = 

Fluorescence  lifetime 

U)p  — 

radius  of  the  pump  beam 

hh.06  = 

radius  of  the  1.06/rm  beam 

^1.32  = 

radius  of  the  1.32/rm  beam 

T*  = 

Transmission  of  the  Q-switch  at  the  lasing  wavelength 

Tg  = 

Transmission  of  the  gain  medium  at  the  lasing  wavelength 

Td  = 

Transmission  of  the  dichroic  at  1.06/iin 

Td  = 

Reflectivity  of  the  dichroic  beamsplitter  at  1.32/xm 

rep_rate  = 

Pulse  repetition  frequency 

Delay  .time  = 

The  delay  time  from  the  time  the  1.32/xm  pulse  is  fired  to  the  1.06/im  pulse 

Pump.pow  = 

Pump  power  incident  onto  the  high  reflector 

V  HR  = 

Transmission  of  the  high  reflector  at  the  pump  wavelength 

VAbs  = 

Absorbtion  efficiency  of  the  laser  crystal  at  the  pump  wavelength 

NpO  = 

Initial  photon  density 

n  — 

population  inversion 

Figure  4.11  provides  a  visual  reference  for  the  definition  of  the  cavity  dimensions. 


All  of  the  length  parameters  are  set  up  in  a  Matlab  program  called  “analyze.m” . 
This  program  passes  the  material,  and  wavelength  information  to  a  program  called 
“qspulse.m”  which  is  the  workhorse  of  the  model.  This  program  calculates  and 


stores  the  photon  number  and  population  inversion  as  functions  of  time.  The  first 
calculation  is  of  the  value  for  ld,  the  distance  the  1.06/zm  beam  travels  in  the  dichroic, 
and  it  is  calculated  from  the  thickness  of  the  dichroic,  t^,  the  index  of  refraction  n^ 
and  the  angle  of  incidence.  The  angle  of  incidence  for  this  optic  was  designed  to  be 
45°.  The  angle  of  refraction  ay,  from  Snell’s  law  is 


o>d  =  arcsin 


sin  (45°) 
nd 


and  from  the  geometry,  the  distance  the  beam  travels  in  the  optic  is 


Id  —  td  COS  (ttd) 


The  next  calculation  is  of  the  effective  cavity  length,  calculated  with  consideration 
for  the  indicies  of  refraction  of  the  various  components.  For  the  1.06/xm  pulse  this 
equation  is 

deff  —  l\  +  riglg  +  l2  +  Tldld  +  h  +  Tiqlq  +  1% 

from  this,  the  amount  of  time  it  takes  a  photon  to  make  one  round  trip  time  is 
calculated 

_  2 deff 

Wt  — 

c 

The  total  loss  in  the  cavity  is  calculated  from  the  transmission  of  the  various  com¬ 
ponents  in  the  cavity 

loss_total  =  (1  —  YlTgTdTqT2) 
and  from  the  above  two,  the  photon  lifetime 

j.  _  trt 

tphoton  —  *  — TT 

loss_total 

The  next  calculation  is  the  Q-switch  rise  time.  The  manufacturer  of  the  Q-switch 
specified  the  turn  on  time  of  the  Q-switch  to  be  “177ns  per  mm  of  beam  diameter 
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inside  the  Q-switch” .  A  function  called  “qsw_rise”  calculates  the  size  of  the  beam 
in  the  Q-switch  from  the  specified  cavity  dimensions.  A  system  matrix  is  calculated 
for  a  round  trip  in  the  cavity  starting  in  the  middle  of  the  Q-switch.  The  spot  size 
in  the  Q-switch  is  then  calculated  using  (35:115) 


L)  = 


\ 


Ml, 2 


7T 


(B.l) 


where  M  is  the  system  matrix.  The  next  calculation  is  of  the  laser  spot  size  in  the 
laser  crystal.  It  is  calculated  using  the  above  equation  except  the  system  matrix  for 
the  round  trip  in  the  cavity  begins  inside  the  laser  crystal.  Using  this  value  of  the 
beam  radius,  the  mode  volume  of  the  laser  in  the  cavity  is  given  as 


V'oo 


TTUl2lg 


2 


The  pumping  rate  gives  the  number  of  photons  absorbed  per  unit  time  in  the  laser 
gain  medium  and  is  given  by 


Pump-rate  =  ump-pow 

he 


The  equivalent  pump  time  is  calculated  from  Eq  (2.24),  repeated  here  for  convenience 


tequiv  —  7/  hi 


where  nj  is  the  population  inversion  left  over  from  the  previous  pulse.  The  initial 
population  inversion  is  then  calculated  from  this  value  and  the  pulse  repetition  rate 
using 

nO  =  Pump_rate  [  1  —  exp  ( teqUiV  - - - 

\  \  rep_rate 
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The  initial  photon  density  is  set  at 


NpO  =  6  x  108 

During  actual  Q-switch  laser  operation  the  power  to  the  Q-switch  is  adjusted  so  that 
the  Q-switch  just  holds  off  the  lasing  action.  This  is  modeled  by  calculating  the  gain 
coefficient  at  the  given  pump  rate  and  calculating  the  necessary  loss  required  to  keep 
the  cavity  from  lasing.  The  gain  coefficient  is 

&SEn  0 
TO  y 
^oo 

the  Q-switch  efficiency  used  to  keep  the  cavity  from  lasing  is  then 

I 

efficiency  =  1  -  -^TJ~  +  0.03 

-L  g-L  d-L  q 

where  the  0.03  is  added  to  put  the  threshold  level  above  the  current  population 
inversion. 

All  of  the  above  calculations  are  necessary  to  set  up  the  calculation  of  the  two 
coupled  nonlinear  differential  equations  which  describe  the  behavior  of  the  photon 
number  and  the  population  inversion.  The  differential  equations  are  given  in  Eq  2.16 
and  Eq  2.17  but  are  repeated  here  (35:253) 


dn 

dT 


-2 — Nphot 

nth 


The  time  in  these  equations  is  normalized  to  the  photon  lifetime.  The  maximum 
number  of  photon  lifetimes  is  set  and  this  value  along  with  the  chosen  integration 
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increment  h  determines  the  total  number  of  loops  that  the  algorithm  makes.  The 
loop  builds  a  picture  of  the  time  changing  photon  density  and  population  inversion, 
one  time  increment  at  a  time.  The  Runge-Kutta  equations  are 


K\  —  /  (n,  Nphot ) 

K2  =  /  +  -h,  Nphot  +  -hKi'j 

K3  =  f  (n  +  -h,  Nphot  +  2^^2) 
K4  =  f  (n  +  h,  Nphot,  hK^) 
the  model  update  is  given  by  the  sum 

RK  —  (-^1  +  2-^2  +  2K3  +  Kf) 


and  it  can  be  interpreted  as  the  average  slope.  These  equations  are  implemented,  in 
Matlab,  in  a  program  called  “runge.m” . 

The  threshold  inversion  level  was  adjusted  for  each  time  step  by  changing  the 
amount  of  Q-switch  loss  according  to 


QjswitchJoss(t)  =  1  — 


efficiency 


This  expression  was  determined  empirically  and  gives  a  90%  -  10%  decay  rate 
equivalent  to  the  value  of  the  Q-switch  switching  time,  r,  used  in  the  equation.  The 
equation  for  the  the  threshold  population  inversion  then  becomes 

Uth{t)  =  2 Igossiy)  ^  (ri  (Q_switchJoss(t)TgT9Td)2  T2)  ^ 

This  expression  is  calculated  in  Matlab  using  a  program  called  “threshld.m” . 
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After  the  program  has  completed  the  specified  number  of  loops  the  pulse  en¬ 
ergy  and  pulse  length  are  calculated  from  the  profiles  of  the  photon  number  and 
population  inversion.  The  pulse  energy  is  calculated  from  the  difference  between  the 
initial  and  final  population  inversions  using  the  modified  version  of  Eq  (2.19)  which 
is  repeated  here  for  convenience. 

Wout  =  ( rii  -nf)  hvL— - 

Qtotal 

where,  again,  the  factor  of  1/2  has  been  dropped  to  account  for  small  lower  level 
lifetime,  n*  and  Uf  are  the  initial  and  final  population  inversions  respectively,  aext 
is  the  loss  associated  with  the  output  coupler,  and  atotai  is  the  total  loss  in  a  round 
trip  of  the  cavity. 
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B.0.1  analyze.m. 


7. 

7. 


analyze.m 


7. 

7. 


'/,  This  program  is  the  master  program  from  which  the  analyses  are  done 

*/. 


clear; 


xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 


xxxx 


These  are  the  INPUT  VARIABLES: 


XXXX 


xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 


7.7.7.'/.  CAVITY  PARAMETERS  7.7.7.*/. 


11=1;  ‘/.cm 

12=4.2;  ‘/.cm 
13=3 . 5 ;  ‘/.cm 
14=7 ;  ‘/.cm 

15=1.7;  ‘/.cm 
16=8 . 6 ;  7cm 


Distance  from  the  input  flat  to  the  laser  gain  medium. 

Distance  from  the  laser  gain  medium  to  the  dichroic  beamsplitter 
Distance  from  the  dichroic  beamsplitter  to  the  1.3  um  Q-switch. 
Distance  from  the  1.3  um  Q-switch  to  the  1.3  um  output  coupler. 
Distance  from  the  dichroic  beamsplitter  to  the  1.0  um  Q-switch. 
Distance  from  the  1.0  um  Q-switch  to  the  1.0  um  output  coupler. 


XXXX 


OPTICS 


xxxx 


Ref_out_10=0.95; 
Ref _out_13=0 . 96 ; 
rad_out_10=20; 
rad_out_ 13=20; 
rad_in=40; 


*/,  Reflectivity  of  the  1.0  micron  output  coupler 
*/,  Reflectivity  of  the  1.3  micron  output  coupler 
*/,cm  Radius  of  curvature  of  the  1.0  micron  output  coupler 
'/.cm  Radius  of  curvature  of  the  1.3  micron  output  coupler 
‘/.cm  Radius  of  curvature  on  the  input  mirror  (999999=>f  lat)  . 


Pump_pow=1.3  ;  '/.W  Pump  power  in  Watts 

Trans_HR=.909;  */,  Transmission  of  the  High  Reflector  at  the  pump  wavelength 
Pump_wav=808 . 5e-7 ; */,cm  Pump  wavelength  in  centimeters 
Pump_spot=0.016;  '/,  cm  Measured  spot  size  (radius)  of  the  pump  beam 


•/.•/.7.*/.  Q-SWITCH  PARAMETERS  XXXX 


pulse_rep_rate=500 ;  '/.Hz 

7,  1.06  micron  pulses. 

Delay_time=200e-6;  '/,<= 

'/,  Q-switch  trigger. 

‘/. 


Pulse  repetition  frequency  of 
Time  delay  between  1.06  micron 

and  the  1.3  micron  Q-switch  trigger. 


•/.'/.•/.*/.  LASER  MATERIAL  WAVELENGTH  AND  POLARIZATION  '/.7,‘/.7. 

material=l;  '/,  Laser  Crystal  material  (l=Nd:YAG,  2=Nd:YLF,  3=SVAP) 

transition=1.0;  l  Lasing  Transition  (1.0  =  1.06  microns,  1.3=1.32  microns) 
polarization=’pi’ ;  7,  Polarization  (only  applies  for  birefringent  host  crystals, 

7,  the  options  are  ’pi’  for  polarization  along 
7,  the  optic  axis  and  ’sigma’  for  polarization 
7.  perpendicular  to  the  optic  axis) 
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Planck’s  constant 
Speed  of  light 


XXXX  PHYSICAL  CONSTANTS  XXXX 
h=6.626e-34;  XJ*s 

c=3e!0;  '/,cm/s 


XXXX  model  fit  parameters  XXXX 

XXX  The  following  parameters  help  the  model  accurately  predict  the  data 


tau_scale_10=l; 

'/,  wo_scale_10=l; 
'/,  Add_loss_10=l ; 
tau_scale_13=l; 

'/,  wo_scale_13=l; 
*/,  Add_loss_13=l  j 


'/,  Scales  the  1.06  micron  q-switch  turn  on  time 
'/,  Scales  the  calculated  1.06  micron  beam  size 
X  adjusts  the  total  1.06  micron  cavity  loss 
*/,  Scales  the  1.32  micron  q-switch  turn  on  time 
'/.  Scales  the  calculated  1.32  micron  beam  size 
'/,  Adjusts  the  total  1.32  micron  cavity  loss 


xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
Xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
xxxxxxxxxxxxxxmxmxxmmmmxxmmm 
•an  calculations  m:i 

xxxxmmmxxxmx 


m:/.  INITIAL  population  INVERSION  XXXX 


n_initial=0; 


'/,  For  the  first  pulse,  it  is  better  to  use  the 
*/.  pulse  repetition 

'/,  rate  variable  to  achieve  the  initial  population. 


•/:/:/:/.  INITIAL  PHOTON  DENSITY  XXXX 


Np0=7e8;  X  photons  WAG 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
xxxxx  xxxxxx 

XXXXX  This  is  the  overall  system  algorithm  XXXXXX 

XXXXX  xxxxxx 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 


total_loss=l; 

XXXX  This  loop  calculates  the  parameters  for  a  series  of  pulses 
for  k=l:10; 

wo_scale_ 10 (k) =l+k* . 1 ; 
for  i=l:16, 
eta_abs(i)=.7; 

pulse_rep_rate ( i) =250+250*i ; 


Add_los8_10  (i) =1 ;  X.1 . 066 ; 

X  Add_loss_13(i)=0.985-i*0.01  ; 
X  wo_scale_13(k)=1.3; 


XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX'/.XX 


xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 


XXXX  For  the  first  pulse  the  pulse  build  up  time  is  set  to 
XXXX  1/2  second  in  order  to  achieve  maximum 
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population  inversion  (this  is  the  last  entry  in  the  call 

7.7.7.X  for  the  function  qspulse. 


7.7.7.  First  pulse  !  ! ! ! ! 
transition=1.0;  7.  microns 
pump_time=l/pulse_rep_rate(i) ; 

[pk.pow, puls_energy,puls_width, prof ile, t.phot ,n_pop,N_th] =qspulse(tau_scale_10, 
wo.scale.10 (k) , total_loss ,11,12,15,16, Ref  _out_ 10 , rad_out_10 , rad_in , Pump_pow , 
eta_abs(i) , Pump _wav, Pump _ spot, material, transition, polarization, n_initial, 1/2) ; 


n_initial=profile(2,size(profile,2))  ;  */,  Leftover  population  inversion  at  the 

*/,  end  of  the  pulse. 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 


7.  transitions. 3; 

7,  pump_time=(l/pulse_rep_rate)-Delay_time; 


7.7.  1.3um  PULSE  1.3um  PULSE  1.3um  PULSE  1.3um  PULSE  1.3um  PULSE  7.7.7.7. 

7.  [pk_pow , puls_energy , puls_width , prof ile , t _phot ,n_pop , N_th] =qspulse (tau_scale_13 , 
wo_scale_13(k) , Add_loss_13(i) ,ll,12,13,14,Ref_out_13,rad_out_13,rad_in,Pump_pow, 
Trans_HR,Pump_wav,Pump_spot .material, transition, polarization, n_initial,pump_time) ; 


7.  n_initial=prof  ile  (2,  size  (prof  ile, 2))  ;  7.  Leftover  population  inversion  at  the 

7.  end  of  the  pulse. 


7<7.7.7.  Now  set  up  for  1.06  micron  pulse 
transitionS.O; 

pump_time=l/pulse_rep_rate(i) ; 

7.7.  1.06  urn  PULSE  1.06  urn  PULSE  1.06  um  PULSE  1.06  um  PULSE  1.06  urn  PULSE  7.7.7. 
[pk_pow , puls_energy , puls_width , profile , t_phot , n_pop , N_th] =qspulse (tau_  scale_ 10 , 
wo_scale_10(k) ,total_loss,ll,12,15,16,Ref _out_10,rad_out_10,rad_in, 

Pump_pow , eta_abs ( i) , Pump_vav , Pump_spot , material , transition .polarization , 
n_initial,pump_time) ; 

n_initial=prof ile (2, size (prof ile, 2)) ;  7.  Leftover  population  inversion 

7.  at  the  end  of  the  pulse. 

7.  transitions .3; 

7.  pump_time= ( l/pulse_rep_rate) -Delay_t ime ; 

7. 

7.7.  1.3um  PULSE  1.3um  PULSE  1.3um  PULSE  1.3um  PULSE  1.3um  PULSE  7.7.7.7. 

7.  7.  [pk_pow ,  puls  .energy ,  puls.width ,  profile ,  t.phot ,  n.pop ,  N.th]  =qspulse  (tau_scale_13 , 
wo_scale_13(k) ,Add_loss_13(i) ,11, 12, 13, 14, Ref _out_13,rad_out_13,rad_in,Pump_pow, 
Trans.HR, Pump.wav, Pump.spot, material, transition, polarization, n.initial, pump .time) ; 

7.  prof  ile_13=  [profiled, :)  ;profile(2, :)  ;profile_13]  ; 

7.  n_initial=profile(2,size(profile,2))  ;  '/,  Leftover  population  inversion 

7.  at  the  end  of  the  pulse 
7.  peak_pow_13(i)=pk_pow; 

7.  energy_13(i)=puls_energy ; 

7.  width_13(i)=puls_width; 
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*/,  phot on_l  if  et  ime_  13  ( i )  =t _phot ; 


Now  set  up  for  1.06  micron  pulse 
transition^.  0; 

pump_time=l/pulse_rep_rate(i)  j 

•/.*/.  1.06  um  PULSE  1.06  um  PULSE  1.06  um  PULSE  1.06  urn  PULSE  1.06  urn  PULSE 
'/.[pk_pow,puls_energy,puls_width,prof  ile.t.phot  ,n_pop,N_th]  =qspulse(tau_scale_10, 
wo_scale_10(k) ,total_loss,ll,12,15,16,Ref_out_10,rad_out_10,rad_in,Pump_pow, 
eta_abs(i) ,Pump_wav,Pump_spot , material, transition, polarization, 
n_initial  ,pump_time) ; 

n_initial=0;  '/,  Reset  leftover  population  inversion 
'/,  for  next  loop 

*/. 

prof ile_10=  [profiled, :)  ;profile(2, :)  ;profile_10]  j  ‘/.Store  the  1.06um  pulse  profile 
*/.  peak_pow_  10  ( i)  =pk_pow ;  '/.  Store  the  1.06  um  peak  power 

energy.  10 d)=puls_energy ;  '/.  Store  the  1.06  um  pulse  energy 

width.  10 d)=puls_width;  */.  Store  the  1.06  um  pulse  width 

photon_lifetime_10(i)=t_phot;  '/.  Store  the  1.06  um  photon  lifetime 

i 

end; 

'/.  Save  data  to  a  file 

data_set=[pulse_rep_rate’ ,width_10’ .energy.lO’ ,  photon.lifetime.10’  ]  ; 
f id=fopen( [’mod.ft’  ,num2str(k) , ’ .dat’] , ’w’) ; 
fprintf  (f  id,  ,a/t1.5e  '/,1.5e  ’/,1.5e  ‘/.l.SeXn’ .data.set’) ; 

fclose(f id) ; 

clear  tau.scale  width.10  energy.10  width_13  energy.13; 
end; 

'/.Save  wo.scale  data  to  a  file 
fid=fopen( [’wo_scale.dat’] , ’w’) ; 
fprintf  (fid,  ’’/.l.Se  \n’  ,wo_scale_10’) ; 
f close(f id) ; 


B.0.2  qspulse.m. 


'/,  qspulse.m  '/, 

7.  7. 

'/,  This  function  returns  the  peak  power,  pulse  energy,  pulse  width,  pulse  profile, 
'/,  photon  lifetime,  final  population  inversion,  and  threshold 
7,  inversion  level  profile. 


function [pk_pow , puls_energy , puls_width , prof ile , t_photon ,n0 , N_th] =qspulse ( 
tau_scale , wo_scale , Add_loss , 11 ,12,13,14, Ref _out , rad_out ,rad_in , Pump_pow, 
Trans_HR , Pump_wav , Pump_spot ,mater_al , trans it_n , polariz , n_ initial , 
pop_build_up_time) ; 


7.7.  mmxxmmmxxmmmxmxmmmm 


7, 7.7.7.  These  are  the  INPUT  VARIABLES:  7.7.7.7. 


7.7.7.7.7.7.7.7.7.7.7.7.7.7.7.7.7.7.7.7.7.7.7.7.7.7.7.7.7.7.7.7.7.7.7.7.7.7.7.7.7.7.7.7.7.7.7. 


'/,11=  Distance  from  the  input  flat  to  the  laser  gain  medium. 

7.12=  Distance  from  the  laser  gain  medium  to  the  dichroic  beamsplitter. 
7.13=  Distance  from  the  dichroic  beamsplitter  to  the  Q-switch. 

7,14=  Distance  from  the  1.3  micron  Q-switch  to  the  output  coupler. 


7.Ref_out=  Reflectivity  of  the  output  coupler 

7.rad_out_10=  Radius  of  curvature  of  the  output  coupler 

7.rad_in=  Radius  of  curvature  on  the  input  mirror  (9999=inf )  . 

7.Pump_pow=  Pump  power  in  Watts 

'/.Pump_wav=  Pump  wavelength  in  nanometers 

7.Pump_spot=  Measured  spot  size  (radius)  of  the  pump  beam 


'/,pop_build_up  time=  time  since  last  Q-switched  pulse 

7.n0=  Initial  population  inversion 

7mater_al  laser  crystal  material 

'/.lambda  lasing  wavelength 

'/.polarization  pi  or  sigma  polarization 

'/.efficiency  Q-switch  efficiency 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

•/.xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx'/. 

7.7.7.7.7.7.7.7.7.7.7.7.7.7.7.7.7.7.7.7.7.7.7.7.7.7.7.7.7.7.7.7.7.7.7.7.7.7.7.7.7,7,7.7,7.7.7.7.7.7.7.7.7, 


7.7.7.'/.  FIXED  VARIABLES  7.7.7.7.7.7. 

7.7.7.'/.  These  are  the  variables  which  are  either  fixed  or  determined  by  the 
7.7.7.'/.  specific  material. 

7.7.7.'/.  OPTICS  7.7.7.'/. 

Ref _in=0. 99983;  7,  Reflectivity  of  the  Input  mirror 
7,  at  both  1.0  and  1.3  microns 

Tg=.999;  7,  One  way  transmission  of  the  gain  medium  (This  turned  out  to  be 

7.  approximately  the  same  for  all  materials,  they  had  good 
7,  coatings  I  think) .  (measured  value) 


B-13 


xxxx 


PHYSICAL  CONSTANTS  XXXX 


h=6.626e-34;  '/,J*s  Planck’s  constant 

c=3el0;  */,cm/s  Speed  of  light 


XXXX  LASER  MATERIALS  XXXX 
XX  XX 


X  Nd:YAG  xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 


if  mater_al==l, 
up_st_lt=237.3e-6;  '/.sec 
Pump_abs=0 . 909 ;  */, 

lg=0 . 5 ;  ‘/.cm 

ng=i.82;  X 

Pump_abs=0 . 9 ;  '/, 

if  transit_n==1.0, 

lambda=l .  064e-4 ;  */,cm 
SE=3.4e-19;  y.cm2 

elseif  transit _n==l. 3, 
lambda=l .  338e-4 ;  y.cm 
SE=0.9e-19;  */,cm2 

else,  stop; 
end; 


Upper  state  lifetime  for  NdrYAG  (measured) 

Absorbtion  of  Nd:YAG  at  808nm  (measured) 

Length  of  the  Nd:YAG  gain  medium 

Index  of  refraction  of  the  gain  medium  at  1.0  microns 
Measured  absorbtion  of  the  pump  power 

Wavelength  of  the  1.0  micron  laser  in  Nd:YAG 
Stimulated  emission  cross  section  at  1.0  microns 

Wavelength  of  the  1.3  micron  laser  in  Nd:YAG 

Stimulated  emission  cross  section  at  1.3  microns 


X  Nd:YLF  xxxxxxxxxxxxxxxxxxxxxxx xxxxxxxxxxxxxxxx 


elseif  mater_al==2, 

lg=0.5;  '/,  cm  Length  of  the  Nd:YLF  gain  medium 

Pump_abs=0 . 99 ;  '/.  Measured  absorbtion  of  the  pump  power 

if  strcmp (polar iz, ’pi’) , 

ng=1.634;  %  Index  of  refraction  of  the  gain  medium  (pi 

*/.  polarization) 

elseif  strcmp (polariz, ’sigma’) , 

ng=1.631;  X  Index  of  refraction  of  the  gain  medium  (sigma 

'/,  polarization) 
else;  stop; 
end; 


if  transit_n==1.0, 

if  strcmp (polariz, ’pi’), 

lambda=1.047e-4;  '/.cm  Wavelength  of  the  1.0  micron  laser  in  Nd:YLF. 

SE=3.7e-19;  '/,cm2  Stimulated  emission  cross  section  at  1.0  microns, 

elseif  strcmp (polariz, ’sigma’) , 

lambda=1.053e-4;  Xcm  Wavelength  of  the  1.0  micron  laser  in  Nd:YLF. 
SE_10=2.6e-19;  '/,cm2  Stimulated  emission  cross  section  at  1.0  microns, 

else;  stop; 
end; 

elseif  transit _n==l. 3, 
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'/, cm  Wavelength  of  the  1.3  micron  laser  in  Nd:YLF. 

'/,cm2  Stimulated  emission  cross  section  at  1.3  microns. 
*/,  This  else  statement  is  here  to  catch  errors. 


lambda=l . 313e-4; 

SE=0.6e-19; 
else,  stop; 
end; 

up_st_lt=503.9e-6;  */,sec  Upper  state  lifetime  for  Nd:YLF 
end; 

W.'l  DICHROIC  BEAMSPLITTER  V/TL 

nd=1.45  ;  '/,  Material  UV  grade  Fused  Silica 

'/,  (Pg  4-13  Melles  Griot  1995-1996  Catalog) 

if  transit _n==1.0; 

td=0. 250*2. 54;  */,cm  Thickness  of  Dichroic 

ang_refr=asin(sin(pi/4)/nd)  ;*/,  Angle  of  refraction  of  beam  into  dichroic 
ld=td/cos(ang_refr) ;  '/,  Distance  the  beam  travels  in  the  dichroic 

T_di=0.929;  '/,  Transmission  of  the  dichroic  at  1.06  microns 

elseif  transit _n==l .3; 
ld=0 . 0 ; 

T_di=.996;  '/,  Reflectivity  of  the  dichroic  at  1.32  microns 

else;  stop; 
end; end; 

mi  q-swiTCHEs  mi 

lq=0.9*2.54;  '/.cm  Thickness  of  Q-switches 

nq=1.69;  '/,  Dense  flint  glass  wavelength  =1.06microns 

prop=.5;  ’/,<=  Used  to  designate  the  distance  into  the 

l  Q-switch  where  the  spot  size  is  to  be  calculted 

'/,  for  the  q-switch  rise  time  calculation. 

if  transit_n==1.0; 

Tq=.975;  '/,  Transmission  of  the  1.0  micron  q-switch 

elseif  transit _n==l. 3; 

Tq=.998  ;  */,  Transmission  of  the  1.3  micron  q-switch 

else;  stop; 
end;  end; 

mi  effective  cavity  length 

l_eff=ll+lg*ng+12+ld*nd+13+lq*nq+14;  %  effective  cavity  length 

TUI  PHOTON  LIFETIME  */.*/. '/,•/, 

t_rt=2*(l_eff )/c;  '/,  round  trip  time 

Add_loss=(l-Trans_HR/4.66)/ (Ref _in*Ref_out*(Tg*T_di*Tq)“2) 
loss_total=(l-Ref_in*Add_loss*Ref _out*(Tg*T_di*Tq)~2)  ’/,  total  losses 

t_photon=t_rt/loss_total;  '/,  photon  lifetime 

llll  q-SWITCH  RISE  TIME  IT/,'/ 
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tau=tau_scale*qsw_rise (wo_scale ,11,12,13,14, rad_out ,rad_in ,transit_n , 

polar iz , prop , lg , ng , lambda) /t_photon ; 


xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 


XXXX  BEAM  WAIST  XXXX 

The  function  cr_spot  returns  the  size  of  the  laser 
beam  waist  in  the  laser  crystal. 


wo=wo_scale*cr_spot (11,12,13,14, rad_out , rad_in , transit _n , 
polar iz, prop, lg,ng, lambda) ; 


xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

•/.xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx'/, 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

XXXX  CALCULATIONS  XXXX 

xxxxxxxxxxxxxxxxxxxx 


PUMPING  RATE 


xxxx 


Pump_rate=Trans_HR*Pump_pow*Pump_wav/ (h*c) ; 

'/,  Assumes  100'/,  transference  to  the  upper 
'/.  laser  level. 

MODE  VOLUME  which  intersects  THE  GAIN  MEDIUM  XXXX 
V00=pi* (wo"2) *lg/2 ; 

XXXX  INITIAL  POPULATION  INVERSION  XXXX 

t_equivalent=-up_st_lt*log( 1- (n.initial/ (Pump_rate*up_st_lt) ) ) ; 

'/,  This  accounts  for  the 
X  leftover  population 
X  after  the  last  pulse. 

nO=Pump_rate*up_st_lt* (1-exp (- (t_equivalent+pop_build_up_time)  /up_st_lt) ) ; 

XXXX  INITIAL  PHOTON  DENSITY  XXXX 

Np0=6.383e8;  X  photons/cm“3  If  this  number  is  too 

X  high  then  the  pulse  starts  to  form 

X  while  the  Q-switch  is  turning  off. 

XXXX  Q-SWITCH  EFFICIENCY  CALCULATION  XXXX 

gamma_0=SE*n0/V00 ; 

gamma_th= ( 1/ (2*lg) ) *log(l/ ( (Ref _in*Add_loss*Ref _out) * (Tg*T_di*Tq) ~2) ) ; 
ef f iciency=l-sqrt (exp (-gamma_0*2*lg) / (Ref _in*Add_loss*Ref _out) ) 

/ (Tg*T_di*Tq) + . 03 ; 
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gamma_th=(l/ (2*lg) )*log(l/( (Ref _in*Add_loss*Ref_out)*( (1-efficiency) 
*Tg*T_di*Tq)“2)) ; 

q-SWITCH  PULSE  CALCULATION  XXXXX 


t_initial=0;'/, 
t_final=75;  '/. 
incr=0 . 1 ;  '/, 

X 


Time  is  in  normalized  units  of  Photon  lifetime. 

Choose  t_final  to  be  after  the  pulse  leaves. 

Size  of  the  differential  step  increment  (too  large  and  the 
whole  thing  becomes  unstable) . 


num_increments=(t_f inal-t_initial)/incr;  X  Total  number  of  differential 

'/,  steps  to  be  calculated. 


prof ile( : , l)=[NpO;nO]  ;  '/,  Row  1  of  profile  contains  the  photon  #. 

X  Row  2  of  profile  contains  the 
'/,  population  inversion 


for  count=2:num_ increments; 

t=  (count-1)  *incr;  I,  Time  increment  in  units 

'/.  normalized  to  photon  lifetime. 
N_th(l,count-l)=(count-l)*incr*t_photon;  */,  build  time  axis 
N_th(2, count-1) =  threshld(t ,Ref _out ,Ref _in,Tg,T_di,Tq,tau,SE,lg,VOO, 

eff iciency,Add_loss) ; 

XX  The  above  line  calculates  the  threshold  as  a  function 
XX  of  time.  It  changes  because  of  the  Q-Switch  loss. 


profiled ,count)=profile(: ,count-l)+  runge(profile( count-1) ,incr,  t,  Ref_out, 

Ref_in,Tg,  T_di,  Tq,tau,SE,lg,VOO,eff iciency, Add_loss) ; 


end; 


profiled,  :)  =  [0:num_iucrements-l]  .*(incr*t_photon)  ;X  build  time  axis 


loss_out=(  1-Ref  _out)  ;  X  Calculate  losses  due  to  output 

7,  coupling 

tUU  Calculate  the  pulse  width  (FWHM)  y.7.7.7.7.7.7.7, 

pk=max (profiled,  :))/2;  7.<=  Find  the  maximum 

puls_width=sum(profile(l, :)>pk)*incr*t_photon;7,<=  Full  Width 


'/.HU  Calculate  the  peak  power  XXXXXXXX 

pk_pow=max(prof ile (1 , :))* (loss_out/loss_total) *h*c/ (lambda*t_rt) ; 

XXXXX  Calculate  the  pulse  energy  XXXXXXXX 

XXXX  The  pulse  energy  is  equal  to  the  change  in  the  upper  state  population 


B-17 


puls_energy= (profile (2 , 1) -prof ile (2 , size (prof ile , 2) ) ) 

*h*c*loss_out/(loss_total*lambda) ; 


return ; 
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B.0.3  runge.m. 


7. 

7. 


7,  runge.m 

7. 

7.  This  function  calculates  an  update  for  one  time  step 

function [RK] =runge(A,h,t ,Ref_out ,Ref _in,Tg,T_dich,T_Q,tau,SE,lg,VOO, 

efficiency ,Add_loss) 

nth_f ull=  threshld (t , Ref _out , Ref _in , Tg , T_dich , T_Q , tau , SE , lg , V00 , 
efficiency, Add_loss) ; 

nth_half =  threshld (t+h/2 , Ref _out , Ref _in , Tg , T_dich , T_Q , tau , SE , lg , 
VOO.eff iciency,Add_loss) ; 

Kl=coup_dif (A,nth_full) ; 

K2=coup_dif (A+(h/2) . *Kl,nth_half ) ; 

K3=coup_dif (A+(h/2) . *K2,nth_half ) ; 

K4=coup_dif (A+h.*K3,nth_full) ; 

RK=(h/6)*(Kl+2*K2+2*K3+K4) ; 

return; 
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B.0.4  coup.dif.m. 


7. 

*/. 

*/. 

% 


coup.dif.m  */, 

7. 

This  function  calculates  the  instantaneous  value  of  the  coupled 
differential  equations  for  a  Q-Switched  laser 


function  [differentials] =coup_dif (A ,nth) 


y immmxxmmxmmmmmmm 

Np=A(l) ; 
n=A(2); 

d_Photons= ( (n/nth) - 1 ) *Np ; 
d_Pop_inv=-n*Np/nth ; 
dif f erentials= [d.Photons ; d_Pop_inv] ; 
return ; 
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B.0.5  crspot.m. 


xxxmxxxmmxxmmxxmmxxxmxxxxxmxx 

'/,  cr_spot.m  '/, 

X  X 

'/,  This  function  calculates  the  spot  size  of  the  laser  beam  in  the  crystal 
function [spot] =cr_spot (11,12,13,14, rad_out , rad_in, transition, 

polar izat ion, prop, lg,ng, lambda) ; 


xxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

XXXX  These  are  the  INPUT  VARIABLES:  '/.'/.'/.'/. 

xxxxxxxxxxxxxxxxxxxx xxxxxxxxx 


'/,11=  Distance  from  the  input  flat  to  the  laser  gain  medium. 

"/,12s  Distance  from  the  laser  gain  medium  to  the  dichroic  beamsplitter. 

'/,13=  Distance  from  the  dichroic  beamsplitter  to  the  Q-switch. 

*/,14=  Distance  from  the  1.3  micron  Q-switch  to  the  output  coupler. 

*/,prop=  Proportion  of  the  distance  into  the 

'/,  laser  crystal  where  the  spot  size  is  desired 

'/,Ref_out=  Reflectivity  of  the  output  coupler 

'/,rad_out_10=  Radius  of  curvature  of  the  output  coupler 

*/,rad_in=  Radius  of  curvature  on  the  input  mirror  (9999=inf ) . 

'/.Pump_pow=  Pump  power  in  Watts 

'/,Pump_wav=  Pump  wavelength  in  nanometers 

*/,Pump_spot=  Measured  spot  size  (radius)  of  the  pump  beam 


'/,Pulse_rep_rate  Pulse  repetition  frequency 

‘/.nO  Initial  population  inversion 

*/,tau  Q-switch  time  constant 

'/.material  laser  crystal  material 

'/.lambda  lasing  wavelength 

’/,lg  length  of  the  gain  medium 

*/,ng  index  of  refraction  of  the  gain  medium 

'/.polarization  pi  or  sigma  polarization 

'/.efficiency  Q-switch  efficiency 


xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 


'/.xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx'/. 


mmmxmxxmmmmmxmxxmxxm 


'/,'/,'/.'/.  DICHROIC  BEAMSPLITTER  %%%% 

nd=1.45;'/.  Material  UV  grade  Fused  Silica 

'/.  (Pg  4-13  Melles  Griot  1995-1996  Catalog) 
if  transition==l .0, 

td=.  250*2. 54;  '/.cm  Thickness  of  Dichroic 

ang_refr=asin(sin(pi/4)/nd) ;  '/,  Angle  of  refraction  of  beam  into  dichroic 

ld=td/cos(ang_refr) ;  '/,  Distance  the  beam  travels  in  the  dichroic 

elseif  transition==1.3, 
ld=0 . 0 ; 
else;  stop; 
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end; 

end; 

•/.m  q-swiTCHEs  xxxx 

lq=0 . 9*2 . 54 ;  */,cm  Thickness  of  Q-switches 

nq=1.69;  '/, Dense  flint  glass  wavelength  ??? 


Tl= [  1  11; 

0  1]; 

T2=[  1  12; 

0  1]; 

T3=[  1  13; 

0  1]; 

T4=[  1  14; 

0  1]; 

T_las_cry_l=[  1  prop*lg/ (ng) ; 

0  1]; 

T_las_cry_2=[  1  (1-prop) *lg/ (ng) ; 

0  1]; 

T_dich=[  1  ld/nd; 

0  1  ]; 

T_Q=[  1  lq/nq; 

0  1  ]; 

Ref _HR=[1  0; 

-2/rad_in  1]  ; 

Ref_out=[l  0; 

-2/rad_out  1]  ; 

M=T_las_cry_l*Tl*Ref _HR*Tl*T_las_cry_l*T_las_cry_2*T2*T_dich*T3*T_Q* 
T4*Ref _out*T4*T_Q*T3*T_dich*T2*T_las_cry_2; 

'/,  The  size  of  the  beam  at  this  point  is  given  by 
spot=sqrt ( (lambda/pi) *M ( 1 , 2) /sqrt (l-((M(l,l)+M(2,2))/2)“2)); 
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B.0.6  qsw-rise.m. 


xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

•/. 


'/,  qsw_rise.m 

'/. 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

X 


'/,  This  function  calculates  the  Q-switch  rise  time  given  the  cavity  parameters, 
function [r ise_t ime] =qsw_r ise (wo.scale ,11,12,13,14, rad_out , rad_in , 
transition .polarization .prop , lg , ng , lambda) ; 


xxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

XX XX  These  are  the  INPUT  VARIABLES:  •/.•/.'/.'/. 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXX 


'/.11=  Distance  from  the  input  flat  to  the  laser  gain  medium. 

'/,12=  Distance  from  the  laser  gain  medium  to  the  dichroic  beamsplitter. 
'/.13=  Distance  from  the  dichroic  beamsplitter  to  the  Q-switch. 

'/,14=  Distance  from  the  1.3  micron  Q-switch  to  the  output  coupler. 
'/,prop=  Proportion  of  the  distance  into  the  laser  crystal  where 
X  the  spot  size  is  desired. 

'/.Ref_out=  Reflectivity  of  the  output  coupler 
'/,rad_out_10=  Radius  of  curvature  of  the  output  coupler 
*/,rad_in=  Radius  of  curvature  on  the  input  mirror  (9999=inf ) . 

'/,Pump_pow=  Pump  power  in  Watts 
'/.Pump_wav=  Pump  wavelength  in  nanometers 

*/,Pump_spot=  Measured  spot  size  (radius)  of  the  pump  beam 


'/.Pulse_rep_rate=  Pulse  repetition  frequency 

‘/.n0=  Initial  population  inversion 

'/,tau=  Q-switch  time  constant 

*/, material  laser  crystal  material 

‘/.lambda  lasing  wavelength 

*/.lg=  length  of  the  gain  medium 

*/.ng  index  of  refraction  of  the  gain  medium 

'/.polarization  pi  or  sigma  polarization 
'/.efficiency  Q-switch  efficiency 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

'/.xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx1/. 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 


nd=l  .45  ;  '/, 

X 

if  transition==1.0, 
td= .  250*2 . 54 ;  '/.cm 


•/.'/.'/.•/.  DICHROIC  BEAMSPLITTER  •/.•/.'/.'/. 

Material  UV  grade  Fused  Silica 
Pg  4-13  Melles  Griot  1995-1996  Catalog 

Thickness  of  Dichroic 


B-23 


ang_ref r=asin(sin(pi/4)/nd) ;  '/,  Angle  of  refraction  of  beam  into  dichroic 
ld=td/cos(ang_refr) ;  '/,  Distance  the  beam  travels  in  the  dichroic 

elseif  transition==1.3, 
ld=0 . 0 ; 
else;  stop; 
end; 
end; 

q-swucHEs  y.ra 

lq=0.9*2.54;'/,cm  Thickness  of  Q-switches 

nq=1.69;  '/.Dense  flint  glass  wavelength 

Tl=  [  1  11; 

0  1]; 

T2= [  1  12; 

0  1]; 

T3=[  1  13; 

0  1]; 

T4= [  1  14; 

0  1]; 

T_las_cry=[  1  lg/(ng); 

0  1]; 

T_Q_2=[  1  ( 1-prop) *lq/ (nq) ; 

0  1]; 

T_Q=[  1  lq/(nq); 

0  1]; 

T_dich=[  1  ld/nd; 

0  1  ]; 

T_Q_1=[  1  prop*lq/nq; 

0  1  ]; 

Ref _HE=[1  0; 

-2/rad_in  1] ; 

Ref_out=[l  0; 

-2/rad_out  1]  ; 

M=T_Q_l*T3*T_dich*T2*T_las_cry  *Tl*Ref _HR*Ti*T_las_cry* 

T2*T_dich*T3*T_Q*T4*Ref _out*T4*T_Q_2; 

'/,  The  size  of  the  beam  at  this  point  (in  centimeters)  is  given  by 
spot=wo_8cale*sqrt ( (lambda/pi) *M(1 ,2) /sqrt (i-((M(l, l)+M(2,2))/2)"2)) ; 
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*/,  The  rise  time  is  given  by 

rise_time=2*spot*177*le-9*10;  */,  See  Manufacturer  specifications 

return ; 
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B.0.7  threshld.m. 


XXXX  threshld.m  XXXX 

*/,  This  function  calculates  the  value  of  the  threshold  inversion 
'/,  from  the  cavity  parameters. 

function [nth] =threshld (t , Eef _out , Ref _in , Tg , T_dich , 

T_Q , tau, SE , lg , V00 , efficiency , Add_loss) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

'/,*/,  t  Time  into  the  pulse  (from  zero) 

'/,*/,  Ref_out  Reflectivity  of  the  output  coupler 

Ref_in  Reflectivity  of  the  input  mirror 

'/,*/,  Tg  One  way  transmission  of  the  gain  medium 

*/.*/.  T_dich  One  way  transmission  (reflection)  of  the  dichroic  045  degrees 

T_Q  One  way  transmission  of  the  Q-switch 

tau  Q-switch  turn  on  time  for  the  particular  beam  size 

SE  Stimulated  emission  cross  section  at  the  lasing  wavelength 

lg  Length  of  the  gain  medium 

V00  Mode  volume  which  intersects  the  lasing  medium 

'/,*/,  efficiency  Q-switch  diffraction  efficiency  at  the  lasing  wavelength 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 


Q_switch=l-ef f iciency/ ( ( (t*l . 6/tau) “3)+l) ; 
gamma_thresh=log(l/(  Add_loss*Ref _out*Ref_in*(  Q_switch* 

Tg*T_dich*T_Q) "2) ) / (2*lg) ; 
nth=gamma_thresh*  VOO/SE; 
return; 
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